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1. Introduction and statement of the theorem 



The nonlinear random Schrddinger equation. 

We seek time quasi-periodic solutions to the nonlinear random Schrodinger equation 

d 

i—u^{eA + V)u + S\u\'^Pu {p>0), (1.1) 

on X [0, oo), where < e, 5 -C 1, A is the discrete Laplacian: 

Aji = 1, \i — j\fi = 1, 

'I -^l^^ (1.2) 

= 0, otherwise, 

V = {vj}j^z'i, the potential, is a family of time independent independently identically 
distributed (i.i.d.) random variables with common distribution g — g{vj)dvj, g e L°°. 
We also assume supp is a bounded set. The probability space is taken to be 

M^'^with measure Y\ giVj) = W 9{'>^j)d'^jj 9 ^ supp ^ is a bounded set. (1.3) 

V = {vjljezd serve as parameters for the nonlinear problem in (1.1). 

Given an initial condition u{0) in £'^{Z^), one of the central questions is whether 
u{t) remains localized for all t, i.e., if tt(0) e £'^{Z^), Vk, can one find R, such that 

\\u{t)\\l2^{Z\[-R,R]y) < Vt? (1.4) 

(Prom now on, we write | | for | || || for || ||^2.) When e = 5 = 0, the answer to (1.4) 
is affirmative. Since u{0) — ^^^^d cijSj, aj — > 0, as \j\ — > oo, u{t) = "^j^id ajdje~^'"i^ 
is almost-periodic (infinite number of frequencies) and the upper bound in (1.4) is 
trivially verified. 

In this paper, for appropriate initial conditions w(0), we construct time quasi- 
periodic solutions to (1.1). So the answer to (1.4) is affirmative for such w(0)'s. This 
is the content of the Theorem and its Corollary. 

Before we enter into the heart of the matter, we first address question (1.4) to 

The linear random Schrddinger equation. 

When 5 = 0, (1.1) reduces to the linear random Schrodinger equation: 

d 

dt ^ ' ' (1.5) 



on X [0, cx)). When < e < 1, it is well known from the works in [AFHS, AM, vDK, 
FMSS, FS, GB, GK, GMP] etc. that the upper bound in (1.4) is satisfied. This is 
customarily called Anderson localization (A.L.) after the physicist P. Anderson [An]. 
Since the potential is time independent: t) = ^(j): properties of time evolution 
can be deduced from the spectral properties of i?, which we summarize below. For 
more details, see the Appendix. 

Let cr(iy) be the spectrum of B.. For E. defined in (1.5), 

cr(iy) = [— 2e(i, 2eci] + supp ^, a.s. (1.6) 

(RecaU the probability space defined in (1.3).) [CFKS, PF]. If < e < 1 and the 
probability measure satisfies (1.3), then almost surely the spectrum of E. is (dense) 
pure point, cr(iy) = cTpp, with exponentially localized eigenfunctions: j £ 

Given u(0) G £^(Z^), we decompose 'u(O) as 'u(O) = ^^^^d Oj^j. So 

u{t) = a,c/.,e— (1.7) 

where Uj are the eigenvalues for the eigenfunctions (^j. u{t) is almost-periodic and 
verifies the upper bound in (1.4). So equation (1.5) has A.L. 

Some motivations for studying equation (1.1). 

Schrodinger equations are equations that describe physical systems, which typically 
correspond to a n-body problem. The linear equation in (1.5) is a 0*^ order app- 
proximation, where the n-body interaction is lumped into the effective potential V . 
Quantum mechanically, |up is interpreted as particle density, so the nonlinear term 
in (1.1) can be interpreted as modelling particle-particle interaction. (The nonlinear 
term in (1.1) can be more general and of convolution type. It will not affect our con- 
struction below.) This is sometimes called the Hartree-Fock approximation (cf. [LL, 
O, Sh]) and is a first order approximation to the original n-body problem. This is our 
first motivation to study (1.1). Other physical motivations along this line appear in 
[FSW]. 

In particular, our method permits us to construct quasi-periodic solutions for the 
Landau-Lifschitz equations on nonlinear classical spin waves with a large random ex- 
ternal magnetic field. Thus 

4- = X [(A5),- + h^] {j e Z'^) 

where Sj are unit vectors in and hj — Vj say; with V — {Vj)j^id a large random 
potential. 
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As explained in [FSW] , we may then seek for a solution Sj ~ 63 and the perturbation 
is subject to an equation of the form (1.1), but with a nearest neighbor convolution 
nonlinearity instead of the local one (see [FSW] for details). As mentioned 

before, (1.1) was chosen as a model but the method described in the paper is sufficiently 
robust to cover in particular any nonlinearity with finite range interactions. 

Our second motivation originates from KAM type of stability questions for infinite 
dimensional dynamical systems. (For results in the standard KAM context, see e. g. 
[E].) (1.1) is a Hamiltonian PDE. It can be recast as the equation of motion corre- 
sponding to a Hamiltonian of a perturbed Z'^-system of coupled harmonic oscillators 
with i.i.d. random frequencies (see (2.2, 2.3)). When 6 — 0, the linear system has pure 
point spectrum: a{H) — a pp. This corresponds to the KAM tori scenario. A natural 
question is the stability of such invariant tori under small (0 < 5 <^ 1) perturbations, 
which leads to construction of quasi-periodic or almost periodic solutions to (1.1). 

Remark. Previously in [AF, AFS], solutions to the nonlinear eigenvalue problem 

(eA + !/)(/) + 5|0|2P0 = ^(/, on ^^(Z^) 
were found, which give the time periodic solutions to (1.1) of the particular form 

u{j,t) = ^{j)e-''^K 

A sketch of the construction. 

We expand in the Fourier basis: e*" '^*(5fe(j) and as an ansatz, seek solutions of the 
form 

«(£,t)= Hj.n)e^''-^%{£), (1.8) 

with the initial condition 

1/ i> 
0) = XI (^kh{£), satisfying ^ \ak\ < 1, (1.9) 

k=l k=l 

where in (1.9), we identify {jk}k^i with {l,...,i/}, Sk with Sj^, {k = 1, ...,iy). The 

unperturbed frequencies are therefore u = cj (V) = V eW, where V = {v^J^i are 
the random potentials at sites jk G Z''. 

Substituting (1.8) into (1.1), we obtain the following equation for the Fourier coef- 
ficients: 

{n-u + eAj + Vj)u{j, n) + 6[{u * vY^ * u\{j, n) = 0, (1.10) 

where v{j,n) — u{j, — n), the convolution * is in the n variable only, *p denotes the 
p-fold convolution and we added the subscript j to operators that originated from 
£^(Z'^). We also write the equation for v: 



{-n ■ LO + eAj + Vj)v{j, n) + d[{u * -D)*^ * v]{j, n) = 0. 
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(1.11) 



Combining (1.10, 1.11), we then have a closed system of equations for y 
we write as 

F{y) = 0. 

Equation (1.12) is a Z''+'^ system of equations. Let j/o = y{t = 0). 

supp I/O = life, -efe}fe=i U {j'fc, efc}fc=i, (1.13) 

where Cfc are the unit vectors of Z^. We seek solutions to (1.12) with y fixed at the 
initial condition on supp yo, i-c, u{jki — Cfc) = a^, t)(jfc5 Cfc) = a^, /c = 1, i/, cf. (2.8). 
We make a Lyapunov- Schmidt decomposition as in [Bl,3, CW1,2]. Let j/o = y{t = 0). 
The equations 

i^(y) =0|zd+.\suppyo on ^^(^d+.'ysuppyo) 

are the so called P-equations, the rest are the Q-equations. The P-equations are used 
to determine y{j,n) on {suppyo}'^- On supp yo, y{j,n) are held fixed at the initial 
condition from (1.9). Instead the v Q-equations determine co = (^(V). 

We use a Newton scheme to solve the P-equations (for more details, see section 3). 
This leads to investigate the invertibility of the linearized operators P/(j/j), where yi 
is the i^^ approximate solution, P/ is P' restricted to [-M*+^, M*+^]'^+'^ (i > 0) for 
appropriate M. 

The random potentials V = {vj^y^^i G are the parameters in the problem. 
Invertibility of Fl{yi) are assured by appropriate incisions in the probability space M^. 
Similar to the linear case in [BW], this is done by using semi-algebraic set techniques 
to control the complexity of the sigular sets and Cartan type of theorem for analytic 
matrix valued functions to control the measure. 

The main difference with the linear case in [BW] is that P/ are evaluated at different 
yi. But due to rapid convergence of the Newton scheme, made possible by estimates 
on Fl,{yi>) for i' < i, this is within the margin of estimates. 

Solving the P-equations iteratively is the main part of the work. The solutions 
to the P-equations are then substituted into the Q-equations to determine u = a;(V) 
iteratively by using the implicit function theorem. Wc obtain time quasi-periodic 
solutions of the form (1.8) to (1.1), which are exponentially localized (both in the 
spatial and Fourier space) to the initial condition (1.9), with modified frequencies 
(jj = (jj{V), which are (e -|- 5)-close to the unperturbed frequencies V = {vjk}k=i- 

We therefore have 



^ ' which 

V 



(1.12) 
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Statement of the Theorem. 



Theorem. Consider the nonlinear random Schrddinger equation 

d 

i—u={eA + V)u + S\u\'^Pu, (p e N+), (1.14) 

where A is the discrete Laplacian defined in (1.2), V = {vj}j^z^ is a family of i.i.d. 
random variables with common distribution g satisfying (1.3). Fix jk G Z'^, k = 
1, • • • ,1^. Let TZ — {jk}k=i C Z'^, V = {vajaen ^ ■ Consider an unperturbed 
solution of (1.14) with e, S — 0, 



k=l 

with lofcl sufficiently small. Let a = {ak}%^i. 

For < e <C 1, 3Xe C R^''\]R'' of positive probability, such that for < S <^ I, if 
we fix X & X^, there exists ^e,<5(x;a) C W , Cantor set of positive measure. There is 
u = u!^^s{V; a) smooth function defined on Ge,6ix', a), such that ifV& Ge,six', a), then 

Ue,sAy,t)= Yl u{j,n)e'^-^'5j{y) (1.15) 
(i,n)ez<i+'' 

is a solution to (1.1 4), satisfying 

u{jk, -Cfc) = ttfc, k = 1, u, 

(i,n)^<S 

|a;- V| < c(e + 5), 

for some c > 0, and where {ek}k^i are the basis vectors for 'Z'^ andS = {jk, —ek}k=i C 
'^d+v _ r/ie sets and Qe,s{x',a) satisfy 

Prob ^ 1, mes R''\ge,s{x; a) ^ as e + 5 ^0. 



Remark. The set C only depends on e; while the set Qe,six;a) C W 

depends on e, S, x E (the random potentials in X^) and a (the initial amplitude). 
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Corollary. There exists X^^s C of positive probability, 0<e<^l,0<(5<^l, 



such that for initial amplitudes a sufficiently small, there are quasi-periodic solutions 
to (1.14). 

Comments on the family of parameters {vj}j^z'i- 

In solving (1.14), we use the basis e*^''^*5j, {j,n) e Z'^+^, (cf. (1.8)). In the Z'^ 
basis Sj (j G Z'^), the hncar operator H = eA + V is not diagonahzed. Hence {vj}jQZ<' 
is not a family of independent parameters. This is a slight variation from the "usual" 
scenario, where the linear operator is diagonahzed and the parameters are independent, 
which is the case in e.g., [B3]. 

Here it is convenient to work with the basis 5j instead of the basis provided by 
the eigenfunctions ipj of H, as i/jj depends on {ffc}fe€Z<i- More precisely, as {vkjkeZ'^XTZ 
is held fixed on the appropriate probability subspace, if^j depends on {vk'}k'eTZj which 
serve as parameters for the construction and are therefore varying, (see the statement 
of the Theorem). 

From the KAM perspective, the normal frequencies are provided by the eigenvalues 
Hj of H. Since {vk}keZ'^\Ti fixed, the strong localization property (A8) of ifjj implies 
that the normal frequencies for \j\ > p, where p only depends on the radius of TZ, can 
in fact be held fixed. This is close to the usual terrain, where the normal frequencies 
are fixed, while the tangential frequencies vary to avoid small divisors, either via the 
parameters or via amplitude-frequency modulation (cf. [B3, KP]). 

An insertion into a larger picture. 

The Theorem presented above is proven for i.i.d. random potentials V — {fjljez^. 
The construction used to prove the theorem is, however general. The essential ingre- 
dient is a spectral separation property on the linearized operator H : H = n ■ uj + H, 
where cv are the tangential frequencies, H is the original linear operator (corresponding 
to the quadratic part of the Hamiltonian, cf. (2.2)). In the present case, H = eA + V. 
Assume H has pure point spectrum and we look at initial conditions localized about 
the origin. Below is a tentative formulation of this spectral property. 

Let pj be the eigenvalues of H. Let i = {j,n), Xi = n ■ uj + pj be the eigenvalues 
of H. Let X be an appropriate function, which depends essentially only on the initial 
condition, localized about the origin, \x\ < 1. Let (f)i, (f)'^ be eigenfunctions of H (i.e., 
products of eigenfunctions of H and the exponentials). Define 



satisfying 



Prob X^^s 1 as e + 5 ^ 0, 
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H has spectral separation property if for each scale L, 3£ <^ L, such that 

|A,-A,H >i^(^,0 (1-17) 

for £ < |i — «'| < L {i ^ i'). The //j, Aj, 0^ can be replaced by their local version 
whenever appropriate. 

In the present case, H = eA + V , we use the local version. Assume e is small so 
that H has A.L. (1.17) is provided by using (A5-7) and restricting to the appropriate 
probability subspace, (2.10) and a direct incision in the frequency space. Related 
spectral separation properties seem to hold in [B3, W]. (Compare (1.17) with the 
nondegeneracy condition in [KP, pi 64], where eigenfunctions do not seem to play an 
explicit role.) 

Remark. For the random Sclirodinger operator H = eA + 1/ (e <^ 1), no Diophantine 
property of the eigenvalues seems to be known at present. So a possible extension 
of the standard KAM method, as outlined in e.g., [FSW] is not feasible. It is known 
following [Mi] however, that the eigenvalue statistics is Poisson and that in a box of size 
N, the eigenvalue spacing is N~'p {p> d). From general considerations, the spectrum 
of H, a{H) is simple [Si]. 

The construction of time quasi-periodic (or almost-periodic) solutions needs a pa- 
rameter. This parameter can sometimes be extracted from amplitude-frequency modu- 
lation, see e.g., [KP]. Nonlinear random Schrodinger equation is an equation endowed 
with a family of parameters, where the separation property (1-17) can be obtained 
from A.L. of the linear operator. So it is a natural candidate for the construction of 
KAM type solutions. 

The continuum Schrodinger equations (linear or nonlinear) are a more frequently 
studied subject. The discrete nonlinear Schrodinger equation presented here should be 
seen as the analogue of the continuum nonlinear Schrodinger equation in a compact 
domain, e.g., on a torus. The lattice therefore can be seen as the indices of the 
eigenvalues or eigenfunctions for the underlying linear Schrodinger operator. 

Time quasi-periodic solutions have been constructed for the continuum nonlinear 
Schrodinger or wave equation in 1-D, on a finite interval with either Dirichlet or pe- 
riodic boundary conditions. See for example, the works of Bourgain, Kuksin, Poschel 
and Wayne in [Bl, KP, Wy]. In [B3], time quasi-periodic solutions are constructed 
for the 2-D nonlinear Schrodinger equation on T^. In arbitrary dimension D, time 
quasi-periodic solutions for nonlinear Schrodinger and wave equations are treated in 
[B5, EK]. 

The construction presented here is related to those in [Bl-5], which use a Newton 
scheme directly on the equations. This direct approach is originated by Craig and 
Wayne in [CW1,2]. It has the advantage of not relying on the underlying Hamiltonian 
structure. The Hamiltonian structure does assure, however, the reality of the frequency 
U3 during the iteration, (see section 2 and [B3]). 
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We end this section by remarking that the present method, as it stands, does not 
yet extend to the construction of almost-periodic sohitions. This is because the hnear 
equation that serves as the starting point of our perturbation is 

i—u = Vu, 
at 

and not 

d 

i^u= (eA + F)M. 
at 

In order to construct almost-periodic solutions, we will need more informations on the 
spectrum of the linear operator ilf = eA + F. 

In [B2] , the construction of almost-periodic solutions for 1-D nonlinear Schrodinger 
and wave equations under Dirichlet boundary conditions was made possible by the pre- 
cise knowledge of the spectrum of the linear operator and the fact that the perturbation 
is quartic (in the Hamiltonian). In the present case it is quadratic. In [B6], almost- 
periodic solutions for a 1-D nonlinear Schrodinger equation under periodic boundary 
conditions and realistic decay conditions were constructed. In particular this applies 
in the real analytic category. Almost-periodic solutions have also been constructed by 
Poschel [P52] in the case of a nonlinear Schrodinger equation, where the nonlinearity 
is "nonlocal". 

PDE's (such as (1.1)) typically correspond to the so called "short range" (but not fi- 
nite range) case. In the "finite range" case, which typically corresponds to perturbation 
of integrable Hamiltonian systems, almost-periodic solutions have been constructed in 
e.g., [CP, FSW, Pol] among others. 



2. Hamiltonian representation and Lyapunov- Schmidt decomposition 

Recall from section 1, the nonlinear random Schrodinger equation 

i^u= {e/^ + V)u + d\u\^Pu (peN+), (2.1) 

where < e, 5 <C 1, A is the discrete Laplacian as defined in (1.2), V — {vj}j^zd 
are i.i.d. random variables with common distribution g as in (1.3). The solutions 

u = {w(j,^)}j6z;'i,te[o,oo)- 

Equation (2.1) can be recast as (infinite dimensional) Hamiltonian equations of 
motion, with canonical variables (u, u) and the Hamiltonian 

=^ Ho{u, u) + 5Hi{u, u). (2.2) 
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Equation (2.1) can then be written as 

iii 



(2.3) 



Remark. The connection with the usual canonical variables (p, q) is u = p -\- iq,u = 
p — iq. The equation of motion in the (p, q) coordinates is 

. dH dH 
oq op 

which can be rewritten as a single equation (2.3). (This also explains the factors i and 
2-) 

Equations (2.2, 2.3) show that (2.1) can be viewed as a perturbed Z'^-system of 
coupled harmonic oscillators with i.i.d. random frequencies. The perturbation Hi can 
be of a more general type, e.g., 

Hi{u,u)= Yl %-^'^r'4^' (2-4) 

with ttjj' = aj'j decaying exponentially or polynomially of sufficiently high degree as 
oo. The reason we mention (2.4) is to stress that the construction we present 
below does not rely on integrability of the system. It also carries through for Hi of 
type (2.4), although we only present it for ajj' = Sjj'. 

The goal of the rest of the paper is to seek time quasi-periodic solutions to (2.1) for 
appropriately chosen localized initial conditions. We hence expand u in the basis 

e'^'-'^'SkU), (2.5) 

where n e Z^,u! e W^, k,j e Z'', 5fc(i) is the canonical basis for Z*^. Skij) is a 
natural basis here due to smallness of e. (In [Bl-3], the spatial basis is given by the 
eigenfunctions of the linear operator. The fc-labeling there is the eigenvalue labeling.) 

In the basis (2.5), (2.1) becomes 

dH 

{n-u + eAj + Vj)u{j, n) + S-^{j, n) = 0, (2.6) 
where n e Z^,j e Z^, Hi as defined in (2.2) and u are the Fourier coefficients of u: 

«(A;,t)= ^«(j,n)e^"-^*5,(A:). (2.7) 
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We have also put the subscript j on operators that operate in the spatial (Z*^) variable 
only. (This is the same notation as in [BW].) 

In view of the Theorem, we seek solutions to (2.6) with the constraint 

w(ife, -efe)=afe {k = l,...u), (2.8) 

where jk € Z'^, are unit vectors in Z'^, are fixed. Assume uji, u)2, ojy are 
rationally independent, i.e., uj = {oJi}^^^ G is a Diophantine vector, which will be 
the case when the Theorem applies. A time shift and a limiting argument (since the 
Kronecker flow is dense) permit us to assume are real. Hence from now on ak G M, 
A; = 1, V. 

Due to the smallness of e, we take our initial unperturbed linear equation to be 

d 

i—u = Vu. (2.9) 
at 

The conditions in (2.8) thus correspond to the initial unperturbed solution 

uo{k,t) = ^a^e-^''^-^%(A;). (2.10) 



Let 

a = {akVk=i e R^ 7^ = {j^jLi C Z^ V - K}ae7^ e K"- 

We constructively show that for e small enough, 3X^ C \R'^ of positive probability, 
satisfying Prob ^ 1, as e ^ 0, such that if we fix a; e X^, for S, a small enough, there 
exists Qe,5ix; a) C R^, Cantor set of positive measure, satisfying mesR'^\^e^5(a;; a) — > 0, 
as e + 5 0. We can find w — uj{V; a), smooth function defined on Qcsi^] ci) cind u 
such that (2.6) holds, uj and u are determined simultaneously in an iterative way. 

Toward that end, we first perform a Lyapunov- Schmidt type decomposition (see 
[Bl-3, CW1,2]) of (2.6). Let 

S = {(j'fc, -ek)\k = 1, . . . z.} c Z'^+r (2.11) 

From (2.10), S = supp uq, uq is a solution to (1-14) when e = S = 0. We call S the 
resonant set and consider the u equations 

dHi 

[{n-u} + eAj + Vj)u]{jk,-ek) + 5-^ {jk,-ek) = ^2.12) 

{k = l,...,iy) 

obtained by taking {j, n) e S. 
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They form the finite system of Q-equations. The remaining infinite system of equa- 
tions are called the P-equations 



dH 

(n • a; + eA,- + V,)^i(j, n) + 5 (j, n) = 0, (j, n) ^ S. (2.13) 

The P-equations are used to determine u{j, n) for (j, n) ^ S. (Recall from (2.8) that 
n), (j, n) G 5} = a are given.) 

Once u{j,n) are determined, the Q-equations in (2.12) are used to determine u = 
uj{V, a) via the implicit function theorem. Since a is real, Hi is a polynomial in u,u 
with real coefficients, the solution u to (2.13) will be real and hence also u determined 
from (2.12). (For more details, see the comment after (2.3) of [B3].) 

To solve (2.13), we duplicate the equation for u to form a closed system. Let 

V = u 
v{j, n) = u{j, -n) 

-S = {{jk, +ek)\k = l,...iy}G Z^+^ (2.14) 

(The flip in sign in the second equation of (2.14) is solely in order that the convolution 
coming from the nonlinearity obeys the usual sign convention.) 

We then have the closed system of P-equations 



(2.15) 



(n • a; + eAj + Vj)u{j, n) + 5 ^ {j, n) = 0, (j, n) ^ S, 
i-n-uj + eAj + Vj)v{j, n) + 5 ^ {j, n) = 0, {j, n) -S. 
For Hi as in (2.2), (2.15) takes the explicit form 

([{n.uj + eA, + V,)u]ij, n) + d[{u * * u]{j, n) = 0, Z^+^\(S U -S) 

\ [{-n-u + eAj + Vj)v]{j,n) + 5[{u*vyP *v]{j,n) = 0, ^"^'"^ ^ 

(2.16) 

where the convolution * is in the n variable only. We solve (2.16) by using a Newton 
iteration scheme to be amplified in the next section. We also identify u with u, u with 

V and write y for 



3. Newton scheme 



Let F denote the left hand side (LHS) of (2.16). Our task is to restrict the set of 
(a;, V) e M^^ in order to find y such that 



F{y) = 0, 
12 



(3.1) 



so that (2.16) is resolved. We use a Newton iteration. Recall first the formal scheme. 

Starting from the initial approximant yo, solution to (1.14) and its conjugate when 
e = S = 0, the successive approximants j/j are defined by 

Ai+^y ''^ y,+^-y, = -[F'{y,)]-^F{y,). (3.2) 
Let T denote the linearized operator F'. From (2.16) 

T = D + SS, (3.3) 
where D is diagonal (in the n & variables) 



D = 



n-(jj + eAj+ Vj 

-n-oj + eAj + Vj 

def f D+ ' 

D_ 



(3.4) 



and 



[p{u*vrp-uv*v {p+i){u*vrp ) > y^-^) 

evaluated along the previous approximant. We note that S is self-adjoint, although 
this does not play a role in our construction. 

Denote by || || the £^ norm of a vector or operator on ^'^{'L'^'^^). Using (3.2), the 
error of the approximation at stage (i + 1) can be estimated 

F(|/i+i) = F(yO + i^'(yO (?/.+!- 2/i) 
+ 0(||y,+i-y,||2) 
= 0(||y,+i-|/,f). (3.6) 

So using (3.2) 

\\ny^+l)\\ = 0{\\[F\y{)]-'f)\\F{y{)\\\ (3.7) 
The crux of the matter is thus to control || [i^'(2/i)]~^ II order that 

\\nyi+i)\\ « wnvM- (3.8) 

(Note the squaring of the norm of F{yi) in the RHS of (3.7), which makes this feasible.) 

Since (3.1) represents an infinite system of equations and the initial condition (2.10) 
is localized in a compact region in Z'*"'"'^, to control the norm of [F'{yi)]~^ we implement 
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the Newton scheme in a shghtly modified way, gradually increasing the size of the 
system that we consider. 

Let M G N"*", which can be assumed large in order that [— M, M]'^'^'^ D 2p suppyo; 
in view of (3.4, 3.5) (See also (3.10) below.). At stage i, let N = M*+^ and TAr(yi), 
the restriction of T{yi) to [-N, N]'^+''. We define 

Ai+iy = y^+l - Vi '^^ - {TNiVi)) ^F{yi). (3.9) 

So 

F{y,+i) = F{y,) + F'{y,){y,+i - y,) + 0{\\y,+i - y,f) 
= {T- TN){y^+l - y^) + 0{\\y,+i - y,f) 

= -[(T - T^)T-']F{y,) + 0{\\T-'f\\F{y,)f), (3.10) 

where we used (3.9). Compared to (3.6) the first term in the RHS of (3.10) is new. 
Moreover it is only linear in F{yi). This necessitates the control of off-diagonal decay 
of T and T^^ evaluated at j/j, in addition to that of ||T^^||. 

The control ofT^'^ 

Recall that the O*'^ approximant (initial condition) to (3.1) yo has compact support, 
suppyo = iS U — iS, where S and —5 are defined in (2.11, 2.14). From (3.3-3.5), T{yo) 
is a diagonal dominated matrix with finite range off-diagonal elements. So off-diagonal 
decay of T{yo) is automatically satisfied. 

Assume the successive approximants yi are uniformly (in i) exponentially localized 
about SU—S (cf. (1.16)). This assumption will be verified later from the construction 
itself in view of (3.9, 3.10). From (3.5) the successive S{yi) have uniformly expo- 
nentially decaying off-diagonal elements in the n direction, and are diagonal in the j 
direction, with a prefactor which decays exponentially in j. (The exponential decay 
of the prefactor stems from the uniform exponential decay assumption on yi.) Hence 
T{yi) are of the type (although more complicated) of matrix operator studied in [BW]. 

To study the T's, we introduce, as in [BW] an auxiliary parameter ^ e R. We 
consider instead 

T^ = D^ + SS, (3.11) 

where 

^ f u-uj + e + eAj + Vj 

\ -{n-oj + d) + eAj + Vj 




(3.12) 



and S is as before in (3.5). Similarly we define Tp^{yi), where = M*+^ as in (3.9). 
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In section 4, fix a; in a good set of probability space, where there is Anderson 
localization for the linear random Schrodinger operator Hj = eA^ + Vj, so that (A6) 
holds. (For precise details see the appendix.) Assuming u Diophantine, yi uniformly 
(in i) exponentially localized about S U —S, we bound the norm of [T^{yi')]~^ , where 
N = M*, i' > i as in (4.9) (The precise relation between i and i' is dictated by the 
construction in section 5.) and establish exponential decay properties of its off-diagonal 
matrix elements on a set of 9 of small complementary measure. 

In Lemma 4.1, fix any y^, we bound [T^(yfc)]~^ for all N. We then use it to 
obtain estimates on [T^(yi')]~^, where = M*, i' > i satisfying the restriction in the 
third line of (4.9). This bound is abstract in the sense that a;, V, yi' are viewed as 
independent parameters for the time being. 

As in [BW], this is an iteration process, using semi-algebraic set techniques and 
Cartan-type theorem. To start the iteration, we neglect 6S and exclude a set of 9 such 
that = has a small diagonal element. To continue the iteration, we also need 
to exclude a set in lv of small measure. It is important to remark that this set in u is 
independent of V, yk- It only depends on x G \E.'^ . We stress that for fixed x in 
the good probability set, Lemma 4.1 holds for any fixed u in the good frequency set, 
any V eM.'^ and any yk which satisfy (Hl-3) in section 4. The set of excluded 9, B{N) 
depends, of course, on x, u, V, yk- 

In section 5, we iteratively transfer the estimates on T^{yi/) in 9 into estimates on 
Tj^{yi) in (a;, V), where N = M*+^, N = M*+^, i > i' > i to he made precise, using 
the resolvent equation and taking into account the Q-equations, which are implicit 
functions relating u, V, yi. (Recall that 9 is an auxiliary variable. In the original 
problem (3.3), 9 is fixed at 0.) x is fixed in the good set of probability space as in 
section 4. 

For the first i iterations, we treat dS as perturbation and use a direct e, 6 per- 
turbation series. Instead of excluding a set of 9 as in section 4, we exclude a set of 
{oj,V) e R^"^, so that in the complement, TN{yi) are invertible with exponentially 
decaying oflF-diagonal elements. 

This generates an initial set of "good" intervals: R^^ D Ai D A2 D • • • D Aj in 

the (cj, V)-space. Using the Newton scheme in (3.9), this also shows that yo, yi, ■ ■ ■ , 
yi are exponentially localized about S U —S. (Recall that yo is the initial condition, 
suppyo = SU -S.) 

Starting from the (z + 1)*^ iteration, aside from direct d perturbation series, for 
certain parts of the estimates, (which concerns the regions far from the origin in the 
direction), we need to keep 5S. This is the heart of the matter. In technical 
terms, we need to deal with more general semi-algebraic sets, which are not solely 
defined by products of monomials. For such semi-algebraic sets, we use Q-equations 
and a decomposition lemma (Lemma 9.9 in [B5], restated as Lemma 5.3 in section 
5) to transfer the measure estimates in 9 in Lemma 4.1 into measure estimates in 
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u) = uj{V) e M.'^ . Using perturbation, this in turn generates a new set of intervals 

Ai_|_i C Ai C ■ ■ ■ C M^^, in the {uj,V) space, on which T^/[z+i{yi) is invertible and 
whose inverse has uniformly (in i) exponentially decaying off-diagonal elements. 

In section 6, using the Newton scheme (3.9), we construct yi-^-i. The (uniformly 
in i) exponential localization about S U —S is preserved. Hence Lemma 4.1 is now 
available at i/k = Vi+i for future iterations. We generate Ai_|_2, yi+2--- 

Section 7 summarizes the entire construction. It is merely meant as a recapitulation 
of the sequence of events leading to the theorem. 

Two technical subtleties 

• The estimates in section 4 are obtained along the construction devised in [BW]. 
However for the application to (4.9), T^{yi>) need to be evaluated at different 
yi' at different scales A^. Due to the uniform exponential decay estimates on j/j/ , 
Lemma 4.1 can be applied as explained after the statement of Lemma 4.1. 

• From the P-equations, the j/j's are constructed on a good set of {uJ,V) G M?^. 
(This set eventually becomes a Cantor set.) On the same set of (cu, V) we also 
have estimates on dyi, where d is with respect to u or V. Using this, we can 
construct yi which is smoothly defined on the whole parameter space (u;,V). 
(Note that outside the good set in (a;, V), yi is no longer close to a solution to 
F{yi) = 0.) Substituting for yi, the Q-equations are therefore defined on the 
whole parameter space (w, V). We can then use the standard implicit function 
theorem to determine uj — uj{V). 



on e^{Z'^+''\{SU-S)), where S, -S are as defined in (2.11, 2.14), which are collectively 



We solve (4.1) using a Newton scheme, with the family of linearized operators T{yi), 
evaluated at the i^^ approximant yi 



4. P-equations and statement in 9 




(4.1) 




T{yi) = D + SS{yi) 



(4.2) 



where 




(4.3) 
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and 

C(„,\^( {p+l){Ui*ViyP p{Ui * Vi)*P-^ * Ui * Ui \ . s 

as in (3.4, 3.5). 

In view of the Newton scheme in (3.9), we need to study the family of restricted 
operators: Tjv(j/i), N = M*+^, M assumed large depending on p, 



TNiVi) ^ RNT{yi)RN (4.5) 

and Rn is the characteristic function of the set [—N, N]'^'^'^. This will be achieved in 
section 5 by using the resolvent identity, covering [— M*+^, _/Vf*+^]'^+'^ with the interval 

/ = [-M\M']'^+'' and smaller intervals J = [-Mo,Mo]'^+^ + /e, |M» < < M*+i, 
Mo ~ (logiV)<^/2, (4.6) 

(see (5.5)), and restricting the set of (a;, V) G M?^. 

Toward that end, as previously mentioned in section 3, we introduce an additional 
parameter ^ e IR and let 

j,e^(D^ + \ 

T\y,) = D0 + 5S{y,), 

= RNT\yi)RN (4.7) 

As mentioned there, we temporarily view a;, V G as independent parameters in this 
section. In the same vein, we also dissociate yi from a;, V, assuming only that they 
satisfy (Hl-3) below. It is only in section 5 that we restrict to a; = uj{V), determined 
from the Q-equations and yj, the i^^ approximate solution to (4.1), which depends on 
a;, V. 

In the rest of the section cu, V are held fixed^ only 9 is varying. Note that is of 

the form '^^T{9) — T'{0 + n ■ to). Later in section 5, we transfer the estimate in 
9 into estimates in a;, hence V by restricting 9 to be of the form 9 = n ■ oj, thereby 
resolving (4.1), (which is at ^ = 0) on the good set of a;, V. 

Newton scheme is an iterative scheme. The estimate on [Ti{yi)]~^ for / defined in 
(4.6) is easily obtained by perturbation arguments on [Tj(yi_i)]~^ known from the 
previous step, which is the step to construct y^, see (3.9). The main task is to estimate 
[Tj{yi)]~^ for J defined in (4.6). We therefore study [T^_j^^ M(|]''+'' (^*-']~^ (and later 
in section 5, we restrict to 9 = k ■ uj^ < 1^1 < M*"*"^). This is the subject of 

Lemma 4. land its application. Note that Mq corresponds to the size of interval at a 
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stage [log Mo/ log M] , which procedes z, while the linearized operator T is evaluated 
aty^: T = T{yi). 

Assume yi satisfies 
(HI) supp Vi C [-M%M*]'^+^ (i > 1). 

(This is by construction, see (3.9).) 
(H2') ||A,2/|| = \\y, - < ^/7T6M-''' {i > 1) 

for some 1 < 6 < 2 in view of (3.6, 3.9) and b will be specified in (6.20). (Recall 

yo = i ]j uo defined in (2.10), vq = uq.) 

\voJ 

(H3) \yi{k) \ < e~"l'^l for some a > (uniform in i). 

(There is no constant in front of the exponential, as we assume small initial data: 
\ae\ < 1,£= 1, ...z/. See (1.9).) 

Remark. Using (3.6) in (3.9) ||Aiy|p < ||Ai_|_iy|| < ||Aiy||, assuming appropriate 
condition on T^^. This is consistent with 1 < 6 < 2 in (H2'). 

(Hl-3) will be verified along the iteration later in sections 5, 6 using Lemma 4.1 below. 

Let An = N]'^, Xn C be a set, where eAj + Vj has A.L. at scale iV, in 

a sense to be made explicit in the process of the proof, Xn is asymptotically (e 0) 
of full measure. (Recall also that V = {vj^.}^^i e with measure Y[k=i9('^jk)^'^jk-i 
is the parameter set.) 

Definition 

For A, c>0, DCa,c{N) C is the set such that if w G DCa,c{N), then 

\\n ■ u;\\j >j^,ne [-N, N]''\{0}. (4.8) 

DCa,c C W is the set of oj such that (4.8) is satisfied for all N. 
Bi3^j{N) C R is the complement of the set of ^ e R defined by 

||T^(y,0-'ll <e^' (0</3<l), 

\[T'M]-\k,k')\<e-^\^-^'^ (7>0), (4.9) 

for \k-k'\> N/10, N = M'+\ i' ~ ^^iKi >i(bas in (H2')), 

log b 

yi' satisfies (Hl-3); i' is chosen in view of later construction in Lemma 5.1, (see in 
particular (5.14, 5.15)). This means of course that at least the first (9( ^°^^ ) approx- 
imants are obtained by direct perturbation series in e, 6. So a = 0(1)| log(e + 5)|. In 
general we write B{N) for Bj3^j{N), unless the parameters (3, 7 need to be emphasized. 
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Remark. At this stage of the construction, it is sufficient to have a lower bound on 
b. This can be easily obtained in the first few perturbation series by adjusting e, 5. 
For later purposes, we mention that the Diophantine condition (4.8) will be used for 
CO — coi', the i' approximation. 

The inequalities in (4.9) are proven iteratively as in [BW]. The rate of decay 7 will 
deteriorate with iteration. So 7 = 7iv- But the decrease will decrease with increasing 
scales and we have 7jv > 7/2 for all N, cf. Lemma 4.1 and the paragraph after. 
This rate of decay determines the rate of decay of j/j. So this is consistent with the 
assumption (H3). 

Inspecting the definition of Dg in (4.7, 4.3), $ is not equivalent to a spectral param- 
eter. Hence we need to resort to Cartan-type theorems as in the wave case in [BW]. 
This necessitates that we obtain estimate (4.9) for more general regions than cubes at 
each scale A^, the elementary regions to be defined below. 

Remark. The various approximants yi are still evaluated using cubes [— A^, A^]'^"'"'^, 
N = M*"*"^ as in (3.9). It is only at each scale that we also look at restrictions of 
T{yi') {i' — ^xo^ i) to these more general regions. 

Elementary Regions 

An elementary region is a set A of the form 

A'^= R\{R + k),ke is arbitrary 

and i? is a hyper-rectangle 

R^{£' e Z^+n|< < Ni,i=l,...d, d + l,...,d+iy, 

£' = {eXi^i e i = {Qi=i e (4.11) 

Let N = maxi Ar/= A^^^^^. Assume £ e fixed. We caU £ the center of R. En{N) 

(at a fixed center) is then defined to be the set of all regions obtained by varying 
k G l/-^'' and A^^ {i ^ imax), keeping A^j < A'". We say 2A'" is the diameter of the 
elementary regions. 

To be economical, we extend the notation T/v to mean TA(Ar) = Ra{n)TRa(n)j for 
any A(A^) G £71{N) and where Ra(^n) is the characteristic function of the set A(A^); 
B{N) is then the corresponding bad set, on which (4.9) are violated. For the purpose of 
constructing approximate solutions, we only need to specialize to N = M*+^, (z > 0). 
However to state the various intermediate technical lemmas, it is more convenient to 
let A^ be any integer. 

Fix any y^ satisfying (Hl-3). Let T/v be the linearized operator evaluated at yk for 
all A^, i.e., T/v = TN^yk)- With a slight abuse of notation, we also let Bp^.y{N) be the 
corresponding bad set. The main goal of this section is to prove 
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(4.10) 



Lemma 4.1. Fix < a < g^qi^jy, a < (3 < 1, Nq sufficiently large, max(^,6(ci + 

v)) <S C < Nq^"^ ■ There exist eo > 0, 5o > 0, such that for aZZ < e < eo, < 5 < 5q, 
there exists X C R^''\R'' with 

mesX > 1 - 0(1)N-'', (4.12) 

where k = k{C,p', d) > 0, p' as in (A2). Fixx e X , 3Q G W (independent ofV G 
Vk), with 

mes Q, < e~^° , 
where k' = k'{C, 13) > 0, such that if 

OJ e DCa,c\^ (4.13) 
then for any A{N) e Sn{N), N > Nq 

mesB^,^^(Ar)<e-^\ (4.14) 
where ^N>Oi- N'^' , (5' >Q), for all N, a as in (H3), a = C>(1)| log(e + 5)\. 

Remark. Bp^^j^i^N) depends only on y^, V as a; is fixed. In the proof of Lemma 4.1, 
only (H3) on j/fc is needed. 

In order to obtain (4.9) at all scales, we apply Lemma 4.1 as follows. From the 
third line of (4.9), for any fixed j/fe, we only need the lemma at scale N = M* with 
i = ^^k. To go to scale N' = = M^^ with the corresponding yy, k' = ^^^C', 
we first use (H2'), which gives that 

\\TUyk)-T'N{yk')\\ < Od+MWvk-yk'W < c»d+.(i)M-^\ (4.i5) 

(4.15) shows that we have essentially the same estimate on ||T^(t//5/)|| as for ||T^(t/fe)||. 

We use as the initial scale instead of Nq, the proof of the inductive step in Lemma 
4.1 then gives instead 

iN'iVk') = MVk) - N-'', {S' > 0), N' = N^. (4.16) 

From (3.9), the decay rate cti+i of yi^i is governed by the decay rate of [T/v(2/i)]~^, 
where N = M*+^. (Note that 9 is fixed at 0.) This operator is treated in section 5 
using several considerations including Lemma 4.1, (4.9), resolvent equation and semi- 
algcbraic sets. The decay rate a^+i depends on 7mo5 where Mq N = M'+^ is 
determined in (5.5), cf. also (6.1-6.8). So (4.16) prevents the deterioration of ai as 
z — > oo and we will have cti > q;/2 for alH in section 6. 
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We prove Lemma 4.1 using iteration. The two pillars of this iteration are semi- 
algebraic set techniques and Cartan-typc theorem for analytic matrix valued functions 
(see [B5, BGS]). The general construction of the iteration is the same as in our previous 
paper [BW]. 

The initial estimate {0^^ step) 

In view of (4.4) and the Newton scheme (3.9), choose M > 2p, such that 

SU-S c[-M,Mf+''. (4.17) 



Lemma 4.2. Fix < a < (3 < 1, there exists Mq — Mo{d + u, a, j3) such that for all 
M > Mq, there exist eo, So such that 

ll[r^,]-^||<e^'' (0</?<l), 

\[TLr\^,^')\ < e-«l^-^'l {a as in {H3)), 
\e-i'\ > M/10,M > Mo, 

for < e < eo, < 5 < So, e e R\Bp^c.{M), m,es Bf3^^{M) < e"^", all x e 
all uj e [0,1)^. (Recall from (4-3-4-5) the x and u) dependence of T^{yl).) 

Proof. We use Neumann series in e, S to estimate T^^. We require 

\ o \ ^ ^ V(i,n)e[-M,Mr^\{5U-5}. 

\~e -n-u; + Vj\ >2max(e-^ ,{e + Sy/^), 

(4.18) 

Clearly, (4.18) holds away from a set of ^ e M of measure at most 

4(2M + 1)^^+^ max(e-^'' , (e + 5)^/^). (4.19) 



Choose e, S such that (e + 5)^/^ < e , which can be satisfied if < e < eo, 
0<S<So with eo = 5o = ^e'^^^ From (4.19), we need 4{2M + lf+''e-^^'' < e'^" 
{0 < a < /3 <1). This leads to M > Moid + z/, a, /3). 

On the complement of the set defined in (4.18), using Neumann series in e,S for 
and (H3), we verify that 

\\[T'm]-'\\ < e""' (0</?<l), 
\[K]-\ej')\<e-^\'-''\ (a as in (ifS)), 

\£-£'\ > M/10,M> Mo, 
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for < e < eo, < 5 < 5o. The probability set at this scale Xm and the frequency set 
at this scale Q^Mi on which and on the complement of which (4.14) holds satisfy 

mes Xm = 1, mes Q,m = (4.20) 

□ 

This direct perturbation argument is the same as in [BW]. Note that (4.20) entails 
that (4.14) holds for all a;,Vj, as in the 0^*^ step, the invertibility is entirely provided 
by shifting in 6. There is no bad site. We will only use the above lemma for the initial 
set of scales. 

The iteration 

We now prove Lemma 4.1 using iteration from scale A^o to Nq = N. [C assumed 
large). We call an elementary region K{Nq) bad, if the first 2 inequalities in (4.9) 
are violated. As in [BW], we need to perform an incision in the frequency space, in 
order that inside any A(A^) e £71{N), there are at most pair- wise disjoint bad 
elementary regions at scale A^o, where means A'"'* (0 < a < 1). For technical 
reasons (cf. [BGS, BW]), this requirement pertains to all elementary regions A(Ar'), 
Nq < N' < 2Nq and not simply at N' = Nq. For later constructions in section 5, it is 
important to note once again that this set is independent of V = {vjk}k=i- 

Let A{N) e £n{N),A{N) be its projection onto Z'^; 

T{N) = {[-A^o, A^o]^ X [-N,N]''}r]A{N), (4.21) 

'T{N) be its projection onto Z'^. Denote by 8TZ{N) the projection of ETZ{N) onto U^. 
Note that A(A'") e ETZ{N) can be of much smaller diameter than 2A'". 

Number of bad elementary regions at scale Nq disjoint from T{N). 

Using (H3) in (4.2-4.4), the region A{N)\T{N) can be treated perturbatively. We 
make separate incisions in the probability space and the frequency space. We first 
make incisions in the probability space. Toward that end, we look at 

A(iV') e Sn{N'), A(iV') c A(A^), A{N') n T{N) = (A^o < AT' < 2No). (4.22) 

Let 

X'^ C M{^\[-^o,iVo]}'* (4.23) 

be the probability set such that there is at most 1 (pair-wise disjoint) A{N') {Nq < 
N' < 2No), satisfying : 

A(A^') eSniN'), 

A{N') C [-N.N]"^, (4.24) 

A{N')n[-No,No]''^(D, 
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where (^1) is violated for some E e I, I = o-{H), the set defined in (1.6). 
Using Theorem A, with m = 7, 



oii){N'^xN;;+^y 

mesXL > ( 1 

No 



> 1 - 



2p' 

0{1) 



^^(.--.(c+D-i) {P'>d{c + 1) + 1), (4.25) 



where 0(1) is a universal geometric constant; p' is to be determined at the conclusion 
of the proof of Lemma 4.1. Wc used N = Nq , the second factor Nq~^^ comes from 
the estimate on the number of elementary regions of sizes 1 to 2A''o associated to 
each lattice site: 0{1) YI'n^'Li -^0'^ (cf- (4.11)) the exponent is an upperbound on the 
number of pairs of elementary regions of sizes up to 2A^o in ^(-^)- 

Given two elementary regions Ai(A^'), A2(A^'), we say that they are convex- disjoint 
if their convex envelops are disjoint. (This is in order that we have (4.23-4.25) at our 
disposal.) To control the number of bad elementary regions at scale A^o, we now make 
additional incisions in the frequency space. Recall that (4.23, 4.24) pertain only to the 
projected elementary regions in Z^. 

We are now ready to prove 

Lemma 4.3. Fix x G X'j^ n Xn^^, where X'-^ is the set defined in (4-23, 4-^4) o,nd 
~ 1 /c 

AjVg is defined as in (4-^3, 4-^4)^ but with Nq' replacing Nq. There exists a set O^, 

mesn^<e-^^, (4.26) 

such that if id <^ O^, then for any A(A) G £1Z{N) any fixed 9, there are at most 2 
convex- disjoint had A{N') G £n{N'), A{N') n T{N) = Nq < N' < 2No in A{N) 
(N — Nq). Moreover is semi-algebraic with degree bounded above by 0{1)N^^'^'^ 
and is contained in the union of at most 0{1)N^^'^'^ components. 

Remark. is independent of V, yk- Observe also that we need localization informa- 

i/c 

tion on the random Schrodinger operators at 2 scales, Nq and Nq' . 

Proof. In view of (4.3, 4.4, 4.7, H3), for A(Ao) such that A(Ao) nr(A) = 0, SS can 
be treated as a small perturbation. We only need to ensure the invertibility of D^^j^^y 
Assume A(Ao), A'(Ao), and A"(Ao) are 3 convex-disjoint 6a(i elementary regions. So 
3{n,j) G A(Ao), {n',j') G A'(Ao), {n",j") G A"(Ao), such that 

\d -\- n ■ uj -\- < 2e~^o or \ - 6 - n ■ uj -\- < 2e~^° ; 

|6' + n' • w + /Xj/ 1 < 2e"^o^ or | - ^ - n' • w + /Xj/ 1 < 2e"^o ; (4.27) 
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and 

\9 + n" ■UJ + I < 2e~^o or \ - - n" ■ oj + nyt \ < 2e~^^ ^ 
where //j, iiy^jijn are eigenvalues of A{No) , A' (Nq) , A."{No) respectively. 
(4.27) implies that 3m, A such that 

|m-a; + A| < 4e-^o , (4.28) 

where m = ±(n — n'), or ± (n' — n") or ± (n — n"), 

X = fij — or — fijff or fj,jff — jjLj. (4.29) 

We use the same argument as in the proof of Lemma 2.3 of [BW], which we summarize 
below. 

There are two possibilities: m = 0, m 7^ 0. When m = 0, from pairwise disjointness 
(4.24, 4.25), (A6) implies 

|A| > e-^o' (0 < < (3) (4.30) 

which contradicts (4.28). 

When m 7^ 0, (4.28) corresponds to at most 

C>(l)Ar^ • (iV^ • N^'^+^f < 0{1)N^'^+'' (4.31) 

number of equations. Since each equation in (4.28) is a monomial of degree 1 in uj, 
the excluded set Q'j^ is of degree less than 0{1)N^^'^'^ . Since |a;| may be assumed to 
be bounded for each such equation, we exclude a set in u of measure O{l)e~^o . It 
is simple to see that for each equation in (4.28), the excluded set in uj has 1 single 
component. 

So in conclusion, for fixed x G X'j^, uj G fi'jy'^' there are at most 2 pair- wise disjoint 
A(A^') G Sn{N'), A{N') nT{N) such that the first inequality in (4.9) is violated 
by using (H3) and a simple perturbative argument. 

Assume A(A^') is such that the first inequality of (4.9) is satisfied. So | ± ^ ± n • 
CO + IJ'jl > 2e~^o for n, j G A(A^') from above considerations. To obtain the second 
inequality we proceed as follows. Let A(A'"') be the projection of A(A^'') onto Z^. In view 
of the restriction in the third expression in (4.9), we may assume diam A(N') > 

We cover A{N') with elementary regions A{Nq'^) of diameter 2Nq'^. 

Since x G X'j^ H Xtvq and on Xtvq for all E (of the form E = 9 + n ■ u here), there is 

at most 1 (pair- wise disjoint) A(A^q' ) in A(A^') where (Al) (with m = 7) is violated. 
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Using the resolvent equation, (Al) and the estimates \ ± 9 ± n ■ uj + iJ,j\ > 2e for 

— 1 /c 

n, j E A{N') for the bad A(iVQ ), we obtain exponential decay in the j direction for 
A{N') for aU E (of the form E = 9 + n ■ oj here). 

We obtain the second inequality in (4.9) for this A{N') by another application of 
resolvent series in 5S and using (H3). This holds for all A(A^') such that the first 
inequality of (4.9) is satisfied, in view of the definition of Xn^. Using A^o = N'^^^ , we 
obtain the lemma. □ 

Number of bad elementary regions at scale Nq intersecting T{N) . 

Wc now estimate the number of bad A(A') G £TZ{N'), A{N') C A(A), A{N') n 
T{N) ^ 0, T{N) as in (4.21), Aq < A' < 2Ao, using semi-algebraic set techniques. 
Here it is important to emphasize the Z'^ coordinate of the center of elementary regions 
as the linearized operator is not a Toeplitz operator in the variable. We look at ele- 
mentary regions with centers in {0} x Z^. We write £TZ{N', j) for the set of elementary 

regions centered at j e Z'^. For any A{N',j) e £n{N',j), let B(N',j) ^= Bp^^iN'.j) 
be a set such that on B'^{N' (4.9) hold. (Later for more general elementary regions 
centered in i e Z'^+^, we will use the same notations as used here.) 



Assume 3Xn' ,ji^N' ,j such that for x e Xn'j, oj e DCA,c{'^N)\Q,N',j, 
mesfi(A',j) < e-^o (Aq < A' < 2Ao, < d < 1). 



(4.32) 



Let 



je[-3iVo,3iVo]<i No<N'<2No 

^No'^= U U ^N',j, 

je[-3No,3No]d No<N'<2No 

^= U U U ^(^''^O' 

j&[-3No,3No]'^ No<N'<2No Sn{N' ,j) 

mes A < 0{l)N^'^+''e-^o from (4.32). (4.33) 

Lemma 4.4. Let N = . For any fixed 9 eR, let 

/={ne [-A,A]'^|n-a; + ^e^}. (4.34) 

Fix X e A^^ . Then for lo e DCA,c{'^N)\n%^, 

\I\ < OdAmo^'^"^ = OaAm'^=N'-'\ (4.35) 
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< 6o < 1 ond assuming 6{d + v) < C < Nq^"^. Hence there are at most 0{1)N^ 
(0{1) a universal geometric constant) pair-wise disjoint bad A{N') e 8TZ{N'), A{N')r\ 
T{N) ^ 0, iVo < iV' < 2No in A{N) (N = N§ ). 

Proof. Since the Green's function is the ratio of two determinants and the norm 
of the Green's function can be replaced by its Hilbert-Schmidt norm, (4.9) can be 
reexpressed as polynomial inequalities in 9. A is therefore a semi- algebraic set. (See 
[Ba, section 7 of BGS].) 

A is defined by 

^ C»(l)7Vo^'^+3^ (4.36) 

number of polynomials, A'^q for number of centers, N^'^^ number of elementary regions 
per center, _/Vq "-^^^^ number of matrix elements. Each polynomial is of degree N^^'^'^^^ 
in 6 (as one squares the matrix elements.) 

Basu's Theorem [Ba], restated as Theorem 7.3 in [BGS], then gives that the number 
of connected components in A does not exceed Nq^'^^^\ If there are n, n', n ^ n', such 
that n, n' belong to the same connected component of A., then from the last inequality 
in (4.33), 

I (n - n') ■ uj\ < C»(l)Ar2'^+^e-^o . (4.37) 

Since n, n' e [-7V, N^^n-n' E [-2N, 2N]'', N = , oj E DCaA'^N) is in contra- 
diction with (4.37) for C < A^q ' < 2 a log Nq (assuming A^o ^ 1), so there can be at 
most 1 integral point in a connected component of A. We therefore obtain (4.35). 

Let i = {j,n) G Z'^+''. Since A{N',i) = A{N',j) + n and T is a Toeplitz operator 
in the variable, we obtain the second conclusion of the lemma. □ 

Remark. C will be a fixed expansion factor. So the upper bound on C will be satisfied 
for all N > No as soon as it is satisfied for an initial A^o- 

A large deviation estimate in 6. 

Lemmas 4.3, 4.4 combined give that the number of bad elementary regions at scale 
ATo in A{N) is at most N'^-^° , where 

N = N^, 6, = l-^(^ [Q{d + v)<C<K/') 

from (4.35). This enables us to use a Cartan-type theorem for analytic matrix valued 
functions (sec [B5]) to prove a large deviation estimate on ||T^(|/j')~-'^ ||, necessary for 
the proof of Lemma 4.1. The proof of the lemma is very similar to the one in [BW] 
(See also [BGS]), after using (4.15) to appropriately adjust j/j' according to the scale 
A^. So we state (without details of the proof) 
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Lemma 4.5. Let bo, (3, a, 7 he fixed positive numbers, so that 

<bo,(3,a <1 and (3 + bo>l + 3a. (4.38) 
Let Nq < Ni be positive integers satisfying 

NoiP, a, 7) < lOOiVo < (4.39) 

with some large constant Nq depending only on P, a, 7. Assume that for any Nq < L < 
Ni, any A(L) e £TZ{L,i), i e Z^+^, 

mes Bf3^^{L,i) < e""^". (4.40) 

Let Xn, be the sets such that on Xn and DCa,c{'^N)\Qn, (4--4-0) holds for all 

1 e [-N,NY+'', L e [Nq,Ni]. LetX'^, Xj^i/c be the sets defined in (4.23, 4.24) and 
Lemma 4-3, the corresponding frequency set as in Lemma 4-3. If 

X e x'^ n Xj^i/c n Xn, 

(4.41) 

a;eL>CA,c(2iV)\{nWuniv}, 

then 

mes{^|||K]-l||>e^'}<e-^'^ (4.42) 

where is restriction ofT to any A e STZ{N), the elementary regions centered at 0, 
provided Nq^ < N < N^'^^ , Ci ^ max(^,6((i+ v)) depending only on P,a. 

Remark. We note that from (4.23, 4.24), the dependence of the probabihty set Xi^i 
(on which (4.40) holds) on z e Z''"'"'^ is only through the Z'^ coordinate. For simplicity, 
we keep the notation i. The set ^ (on the complement of which (4.40) holds), 
on the other hand, does have full dependence on i. 

A summary of the proof. 

We use analytic and harmonic function theory together with a 2-scale (in the range 
[A^o,-^i]) analysis to control the measure of the set in (4.42) at scale N S> A^o- (See 
the proofs of Lemma 4.4 in [BGS] and Lemma 6.2 in [BW].) For a given N, let these 

2 scales Li, L2 G [A'o, A^i], Li < L2, satisfy 

log Li ~ log A^, log L2 ~ - log Li ~ — ^ log A", 

C\ as in the lemma, Ci ^ 1 /u. We reiterate the main line of arguments below. 
• Fix 6*. Let 

A^ = {m e A(Ar)|3Ai e £R{Li),Nq < Li < 2Aro,Ai c m + [-Li, Li]'^+'^, Ai is bad}, 

(4.43) 
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For X, u) satisfying (4.41), 



\K\<N^~^° (60 >0), (4.44) 

by Lemmas 4.3, 4.4. Since X'^^ D X^, fl'^^ C O7V, where X'^^, Q'^^ as defined 
by the first two equations of (4.33). 

• Let A* be, roughly speaking, the complement of the set in (4.43). For more 
precise definition, which requires a partition of A, see the beginning of the proof 
of Lemma 4.4 in [BGS]. Using an elementary resolvent expansion. Lemma B in the 
appendix, which is Lemma 2.2 of [BGS] reiterated, we obtain an upper bound on 
ll[^A»]~^ll fixed 9 by using the decay estimate on the Ai's, elementary regions 
at scale Li, in A*. By definition they are all good. By standard Neumann series 
arguments, this bound is preserved inside the disk B{9, e~^°) C C. 

Remark. We have control over the size of A* via (4.43, 4.44), but not its geometry. 
Typically A* is non-convex, hence the need for elementary regions which are more 
general than cubes, in particular L-shaped regions, in view of Lemma B. 

• Define a matrix- valued analytic function A{0') on B{0, e~^°) as 



D rji9' p [1^0' 1 ^ D T'^' T? 

- -n-As-tjv-n-Aj-'Aj -n-A^-t AT-n-As 



(4.45) 



where A^ = A\A*, -Ra,, -Raj are projections. From (4.44), A{0') is a rank 
C)(iV^~^") X 0(iV^~''o) matrix. The raison d'etre of introducing A{9') is the 
following inequality: 

P(^')-'ll<ll(<)-'ll<e'^°P(^')-'ll, 

(see Lemma 4.8 of [BGS]). So to bound IK?"^)"^!!, it is sufficient to bound 
||A(0')~^||, which is of smaller dimension. 

Toward that end, we introduce an intermediate scale L2, logL2 ~ (logLi)/a" > 
logLi. We work in an interval O ^ {9'\\9' -9\ < e~^"}. Using (4.40) for the A2's 
at scale L2 and in A, Lemma B, we obtain an upper bound on [[[T"^]""*^!! except 
for a set of 9' of measure smaller than e"^^^^). So there exists y G 0, such that 
we have both an upper bound on ||yl~-'^(^')|| at 9' = y, hence a lower bound on 
the smallest eigenvalue of A{9) and an a priori upper bound on ||yl(^')||, which 
comes from the boundedness of and the bound on [T^*]"^ (see (4.45)). 

Transfering the estimates on ||y4(^')||, ||A~^(^')|| into estimates on log | detA(^')|, 
which is subharmonic and using either Cartan-type theorem (see sect. 11.2 in 
[Le]) or proceeding as in the proof of Lemma 4.4 of [BGS] or Chap XIV of [B5], 
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we obtain the lemma by covering the interval / = {-0{N['^^),0{N^^^)) with 
intervals of size e~^°. (Recall (4.39) and that for all 0^1, is automatically 
invertible.) □ 

Iteration Lemma 

To obtain exponential decay of T^^, we need 

Lemma 4.6. Suppose M, N e N"*" are such that for some < r < 1 

<M < 2N^. (4.46) 

Let Aq G £1Z{N) he an elementary region with the property that for aZ/ A C Aq, A e 
8n{L) with M <L<N, 

\\{Tl)-^\\<e^' (0</3<l). (4.47) 

We say that A e £TZ{L), A C Aq is good, if in addition to (4-4V} 

\\{Tl)-\k,k')\\ <e-^l'=-'='l (4.48) 

for all k' E A,\k — k'\ > L/10. Otherwise A is called bad. 

Assume that for any family T of pairwise disjoint bad M' -regions in Aq with 

M+1<M'<2M+1 

< N^{0 < l3 < 1). (4.49) 

Then one has 

|(T^)-^(A;,A;')I <e-^'""''' (4-50) 
for all k, k' e Aq, \k - k'\ > 7V/10, 7' = 7 - N'^' {5' > 0), provided N > 7Vo(/?, t,7). 

The proof of the above lemma is written out in detail in [BGS]. The only difference is 
that instead of being tridiagonal, has exponentially decaying off-diagonal elements 
from (H3). So 7 = 7(0;). The proof goes through. So we do not repeat it here. (See 
also [BW].) The gist is as follows. 

The exponential decay estimate at scale N in (4.50) is obtained from exponential 
decay estimate in (4.48) at smaller scales M' by using (4.15), the norm estimate in 
(4.47) and the resolvent identity. To implement this, we use a sequence of scales 
Afj+i = Mf\ with Mo = M, C is chosen such that C 13 < 1, C'r < 1. For each 
elementary region at scale Mj_|_i, A(Mj_|_i) and for each k G A(Mj_|_i), we exhaust 
A(Mj_|_i) by an increasing sequence of annuli centered at k of width 2Mj, or more 

29 



-N^ (4.51) 



precisely the intersection of this sequence with A(Mj+i). Roughly speaking, an annulus 
is good, if it does not intersect a bad cube of the previous scale Mj. 

In each of the connected components of the complement of the bad annuli, we apply 
the resolvent identity using the estimate in (4.48) for elementary regions of size Mj. 
In the bad annuli, we use (4.47). From (4.49), the number of bad annuli is at most 
sublinear in Mj^i. Using a multiscale induction argument to reach the scale AT, we 
obtain the exponential decay in (4.50), when \k — k'\ > N/10. 5' is determined from 
(4.47, 4.49),5'~t(1-/?C"). □ 

Let Bi3^j{N,i), i e Z'^+'^ be a set such that on the complement, (4.9) hold. When 
z = 0, we write Bf3^-y{N, 0) = Bi3^^{N). As before let 

XN,i, ^N,i be the probability, frequency subsets on which and 

on the complement of which mes B/3,^(A/', z) < e 

Combining Lemmas 4.5, 4.6 with (4.23-4.25), Lemmas 4.3, 4.4, we obtain 

Lemma 4.7. Assume that for any No < Nq < Nq , max(^, 6(ci + v) <^ C < Nq^'^ 
(Noie.S) > 1), any A{No) E £n{NQ,i), i e 2^^+-^ 

mes-B^,^(Aro, i) < e'^^ {0 < a < 1). (4.52) 

Let [Nq, N^^''] be the next interval of scales. For any N e [N^ , N^''''] , write N = N^ 
with Nq e [No,Nq]. Let X'j^, Xj^i/c be the sets defined in (4-23, 4-^4) osnd Lemma 
4-3, Q.'pj the set defined in Lemma 4-3, satisfying 

-^^^^ ^ 1 - j,Mp^-iiL)-i) > + 1) + 1)' 

mesX^vc > 1 - iP > d{C- + 1) + 1), (4-53) 

JL 

ml ^ —N2C 

mesU^ < e , 

in view of (4-25, 4-26). 

Let Xn, be the sets such that on Xn and DCa,c{'^N)\Qn, (4-52) holds for all 
i e [-Nf%Nf-Y+-, all No e [No,N^]. If 

X G X'j^ n X^i/c n Xn = Xni 

mes Xn > 1 — with p" = p"{p', C, d + u) > 1 by choosing p' large enough; 

CO e L>CA,c(2Ar)\{0'^ U Aat} =^ DCA,cm\nN, (4.54) 

mes iiN < e , 

Q,N is semi- algebraic with degree less than A"*^ (d+i^) 
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then 

mesB^,y(7V) < e"^" (0 < (t < 1). (4.55) 
where 7' = 7 - N'^' (5' > Oj. 

Proof. Applying Lemma 4.5 using (4.52), we obtain the large deviation estimate on 
||[T^]~^||. Choosing < r < for x, uj in the sets defined in (4.54), (4.49) is 
satisfied. (Here we need the definition that an A'^^-region is bad if it intersects a 
bad A^o-region. Otherwise it is good. On the good A^'^-region, (4.48) is obtained by- 
using resolvent expansion (Lemma B) and (4.9) for A'o-regions.) Hence Lemma 4.6 is 
available and we obtain (4.55). The estimates on Xjv, follow from (4.51, 4.53) and 
the constructions in Lemme 4.3-4.5. □ 

Proof of Lemma 4.1. Assume Nq (to be determined below) is such that Lemma 
4.7 is available. For the scales Nq < N < Nq , we use Neumann series in e and 5 a la 
Lemma 4.2 and its proof. For the scales A" > A^g^, we use Lemma 4.7. 

From Lemma 4.2, we need 

Nq> MQ{d + u,a,(3) (4.56) 

and 

eo = SQ = le-^^o\ (4.57) 
From Lemma 4.7 and the choice of cr, the expansion factor C needs to satisfy 

max(-,6(d + zy)) < C < iVo^/^ (4.58) 
cr 

From Theorem A and (4.19), Nq needs further to verify 

iVo > max((5,Mo), (4.59) 

where Mq = Mo((i + u,(3,a) as in Lemma 4.2. Fix A^o satisfying (4.59), C satisfying 
(4.58). (4.57) determines eg, Sq. 

For the scales Nq < N < A"^, A^o satisfying (4.59), the estimates 

||[T^]-^||<e^' (0</?<l), 

\[Tn]~\^J')\ < e-"l^-^'l {a as in (Jf3)), 
\e-e'\ > N/10, 

for < e < eo, < (5 < (^o- ^ire obtained using Neumann series by shifting in 9 only. 
So 

mesXAr,i = 1, mesnN,i ^ 0,{Nq < N < N^), 
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where Xjv,i, ^N,i are as defined in (4.51), and (4.52) hold. 
Let 

N je[-3Arc,3Arc]<* 

where the second equahty follows from the remark after Lemma 4.5, j being the 
coordinate of i; 

^=U U (4-61) 

N ie[-3JVC,3JVC]'*+"' 

where Xnj, flN,i are as defined in (4.51). 

On X, DCa,c\^j Lemma 4.7 is available with the initial 7 = 0: from Lemma 4.2 
for iteration to all scales. Estimating the measure of X and O using (4.53, 4.54, 4.60, 
4.61), using the measure estimates on the bad set in 9 from Lemma 4.2 and 4.7, we 
obtain the Lemma by taking p' ~ 0(i{l)C^- D 



5. Invertibility of T{yi), Q-equations and determination of u 

Fix X E X C ]R^''\]R'^, defined in (4.12), which generates a corresponding set O as in 
Lemma 4.1. The main work of this section is to convert the measure estimates in 6 for 
fixed uj G R'^\0, fixed V G M'^ in Lemma 4.1 into measure estimates in a; = uj{V) G M.'^ 
and extend to (w, V) G M^'^ in the neighborhood of a; = i^(V), while keeping 9 fixed: 
^ = and addressing the original family of linearized operators Tjv(yi) defined in 
(4.2-4.5), where yi is the i^^ approximant to the P-equations in (4.1). This is possible 
because is only a function of n • a; + ^, = T'{n ■uj + 9),cL (4.2-4.4, 4.7). 

Since 9 is now fixed at 0, before doing the conversion, we need to make a further re- 
striction in X in order that the spectrum of the various restricted random Schrodinger 
operators stay away from 0. This is needed when n — and we cannot vary u to have 
invertibility of the linearized operators, cf. (4.3). 

So we modify the definition of X'j^ in (4.23, 4.24) to also include the condition 

dist {a{A{N')), 0) > e-^o' (0 < < /3) (5.0) 

for aU A(iV') n [-N,N]'^ ^ 0, {Nq < N' < 2iVo), i-e., we require (4.30) to hold also 
when A = /Uj, eigenvalues of restricted random Schrodinger operators and not just 
differences of pairs of eigenvalues. In view of the Wegner estimate (A7), this leaves 
the measure of the set X'j^ in (4.25) essentially unchanged. 
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This generates the restricted probabihty set X C X C \]R'^, on which Lemma 
4.1 holds. Rename X as X and let 0, be its corresponding frequency set. 

To do the conversion, we need to supplement the measure estimates in (4.14) by the 
fact that the bad sets B{N) defined from (4.9) has a semi-algebraic description in terms 
of (w, V, 6*), enabling us to use the decomposition Lemma 9.9 of [B5]. Once we have 
the necessary estimates on [Tjv(j/i)]~^ after removing a small set in = i^i(V) (we 
put the suffix i here to stress that it is the i^^ approximation), we construct the next 
approximant j/j+i according to (3.9), which in turn is used to construct uJi+i = uJi+i (V). 

In this section, we primarily address the invertibility of Tis[{yi), N = M*"''^. Since 
the estimate on TN{yi) and the construction of j/j+i are interconnected, it is good 
at this point to lay down the complete induction hypothesis. (In section 4, we used 
first part of the induction hypothesis (HI, 3) and the first inequality of (H2) to derive 
Lemma 4.1.) The first few approximations are obtained by using direct perturbation 
series in e, 5, (see (4.9), as well as the text and the remark afterwards). So < e < 
eo < 1, < (5 < (5o < 1, and we have a = 0{\)\ log(e + 5)\. 

On the entire (a;, V) parameter space, we assume: 

(HI) suppy, C [-M\M^Y+'' {i > 1), 

(H2) ||A,y|| < \\dAiy\\ < 5„ 

where || || stands for sup^ || (Recall that we identify y with y; d refers to 

derivation in a; or V.). 5i,5i will be shown to satisfy: 

5i < \/e + 5M"^^3)\ 
(H3) \yi{k)\ < e-"l'=l for some a > 

(The constant in front of the exponential in (H3) is 1, because \a(\ <^ 1, £ = 1, 
see (1.9).) 

From (H2), y is aC^ function on (oj, V). Application of the implicit function theorem 
to the Q-equations in (2.12) with 




yields 

= V + (e + ^)^/7^(V), with Wd^^W < C, (5.2) 

whose graph we denote by F^. Recall that a priori, yi are only defined on certain 
intervals in {u, V) space. It is in order to use the implicit function theorem that we 
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extend yi to the entire {ciJ,V) space, using the estimates on dyi in (H2). (H2) and 
(2.12) imply, moreover, that 

IIV'^ - V'^-ill < 0{l)\\yi-yi-i\\ < 0{l)6i, (5.3) 

which in turn imphes that 

Fj is an (e + 5)5j-approximation of rj_i. (5-4) 



At each stage i, define 

Mo = 0{l){i + if/'' {log Mf/^ (5.5) 

for some C > and ^Mo,k as defined in (4.51). Unhke (Hl-3), the following hypothesis 
is only assumed to hold on certain intervals in M?^, the (a;, V)-parameter space. 

(H4) There is a collection Sj of intervals / in M?" of size M"**^, same C as in (5.5), 
such that 

(H4, 0) (inVi) C {DC'^l\Qi), where 

DCf/= DCA,ci'2N), N = M\ 

def 

= U;jezdn[-2Mo,2Mo]'* ^Mo,k, i > iQ > 

(see the remark after (H4, iv) concerning iq). 

(H4, i) On / e T,i,yi{(jj,V) is given by a rational function in (a;, V) of degree at most 
Mi^\q e N+). 

(H4, ii) For {u,V)e\Jie^J 

\\dF{yi)\\<Ri. (5.6) 



In (6.20) Ki^ Ki will be shown to satisfy: 



K 



J 1/4 \i+2 



(5.7) 



(H4, iii) For (a;,V)eU/eE,^, 

TN^TN{y^-l),N^M' 
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satisfies 



\T-\k,k')\ < e-^l'^-'^'l for \k-k'\ > 
(H4, iv) Each / e Ej is contained in an interval /' e and 

mes.(r,f|( U A U n)<M-'/'. 
/eEi_i /eEi 

Remark. (H4, 0) is only needed for i > {q to ensure the availability of Lemma 
4.1. Up to stage io, we use direct e, S perturbation series, where the Diophantine 
property of u> is not required, (cf. Lemma 4.2). 

Contrary to related estimates on Ajy,5Ajj/ in (H2). (5.6, 5.7) cannot be extended 
to the entire {co, V) space. This is because, as mentioned earlier, outside the intervals 
in Sj, Hi are no longer close to solutions to the P-equations in (2.16). 

Invertibility ofTNiVi), N = M^+^ 

Assume (Hl-4) hold at stage i. To construct j/j+i, we need to control 

[TN{y,)]-\N^M^+\ 

with a further restriction of the (a;, V)-parameter set. This will give us (H4,iii) at stage 

Wc accompUsh this by covering [— M*+''^, *+^]'^+'^ with [— M% M*]'^+'^ and intervals 
[-Afo,Mo]^+^ + /c. Mo as in (5.5), k E Z^+^, MV2 < \k\ < M^+^ and using the 
resolvent identity. We first estimate [TMi{yi)]~^ ■ Fix (o), V) G U/es (H4,iii) at 
stage i gives 

||[TM.(y^-i)]-i <M^'', 
\[TM^{y^-l)]-'{k, k')\ < e-«l^-^'l (|A: - k'\ > i^). 



(5.8) 



We write 

TMiiVi) = TM^iVi-i) + [TMiiVi] - TMiiVi-i)] 



(5.9) 



Prom the first inequality in (5.8) 



|A-^||<M*''. (5.10) 



The first inequality of (H2) at stage i gives 



\B\\<0{l)M-^^^\ (5.11) 
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So 

\\[TMi{yi)]~'\\ < 2M^^ for i > C\ (5.12) 

To obtain pointwise estimate on T^](yi), we use (5.9) and resolvent series. has 
off-diagonal decay from (5.8), B has off-diagonal decay from (H3) at stage i. Iterating 
the resolvent series and using (5.12), we obtain 

\[TMi{yi)]-\k,k')\<e-"'\''-'^'\{\k-k'\>i''), (5.13) 

with a' = a — M~'^ > a/2 (d' > 0), uniformly in i. 

We now study |[T(yi)]~^| on the Mq intervals, Mq as in (5.5). We distinguish 2 
types of Mo intervals J in [-M^+\ M^+^+^r 

• J n [-Mo, Mo]*^ X [-M'+\ M'+^]'' = (a) 

• Jn [-Mo,Mo]'^ X [-M*+i,M»+i]^ 7^ (6) 

For type a, we use direct perturbation in view of (H3). For type b, we use a more 
delicate construction. We write Mq = M*°. (Mo is chosen in order that the total 
degree of the semi-algebraic set describing the bad set in a;, V, 6* is not too large.) 

The Mq intervals are at [ioY'^ scale. On Mo intervals of type 6, we use Lemma 
4.1, which is for TMoiyio)^ — ^^y^o < ^ similar to (4.9). Using a decomposition 
lemma ([Lemma 9.9, B5] restated here as Lemma 5.3) to make appropriate incisions 
in the {u>, V)-parameter space, applying (H2) between j/j^ and i/i and combining with 
estimates on type a intervals, we obtain 

\\{RjT{y,)Rj)-'\\ < e^o (0 < /? < 1) 

M (5-14) 

\{RjT{y,)Rj)-\k, k')\ < e-«"l^-^'U, k' eJ,\k- k'\ > ^, 

where a" ^ a - M''^" {6" > 0), for aU 

J = [-Mq, Mo]'^+'' + k, \m' < \k\ < M'+\ (5.15) 

This is the content of Lemma 5.2. We delay its precise statement and proof momen- 
tarily. We first prove 

Lemma 5.1. Assume (5.12-5.15) hold. Mq is as in (5.5). Then 

\\[TMi^.{yi)\-'\\ < 0{1)M^'^^)'' (5.16) 

\[TMMy^)r\k,k')\<e-^\^-^'\ for\k-k'\ > {i + if , (5.17) 
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with a = a- M-(*+i)^, ^ > 0. 

Proof. (5.16, 5.17) are exercises in resolvent identity or equivalently using Lemma B 
in the appendix. We first prove (5.16). (5.17) then follows by using (5.16) and another 
application of the resolvent identity. Let 

For any k,£ E -Bj+i, we have (Assume T^^^^ is defined.) 

= Tw\k)if'^^)+ E T^\k)ik^f'')TBl,{k",e) (5.18) 

k'ed^W{k) 
\k"-k'\ = l 
k"€Bi+i\W{k) 

for \k\ < ^M\W{k) = Bi, for \k\ > \M\ W{k) is a size Mq interval. It is easy to 
see that V/c', 3W{k), such that disi{k' , d^W {k)) > Mq, where d^,W{k) is the interior 
boundary of W^(A;), relative to Bi+i. Summing over £ e -Bi+i yields 

,s^P E \TBUk,i)\< sup Yl \\^w\k)\\ 

+ sup Yl \Tw\k)ik,k')\ sup Yl (5-19) 

Using (5.12-5.15), we have 

sup Y l^eii ^) I ^ ^M*'' • (2M* + 1)^^+" 

+ 0(l)e-=^°Mo''+^-^ sup Y \TB'^,ik",e)\, (5.20) 



since > Mq. So 



I -Di+l II — 



which is (5.16). 

To obtain (5.17), we retrace our steps back to (5.18) and restrict to k,£, such that 
\k — i\ > {i + 1)*^ ^ Mq, in view of (5.5). Iterating (5.18) along the path from k to £ 
using Bi, J and using (5.16) for the last factor, wc obtain (5.17). □ 

(5.16, 5.17) will be the conclusion of (H4, iii) at stage i + 1 once we specify the 
new set of intervals on which they hold. As alluded to earlier, will be 

determined from and the new restriction on (w, V) in order that (5.14, 5.15) hold. 

Determination of Sj+i 
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Lemma 5.2. Assume (Hl-4) at stage i. There exists f j C Fj, meSj^f j < M such 
that (5.14, 5.15) hold on 

U (/n(rAf,). (5.21) 

The proof of (5.21) relies on the measure estimates in Lemma 4.1 and semi-algebraic 
description of the bad set. We need the following decomposition lemma, which is 
proven in [B5, Lemma 9.9]. 

Lemma 5.3. Let S C [0, 1]^"' he a semi- algebraic set of degree B and mes2n»5 < 
r7,logi? -C logl/ry. Denote by {x,y) G [0,1]" x [0,1]"^ the product variable. Fixe > 
^i/2n j'/jgj^ there is a decomposition 

with S\ satisfying 

\Proj^Si\ < e (K > 0) (5.22) 
and S2 satisfying the transversality property 

meSn(cS2 n L) < B^e'^rj^/^'' {K > 0), (5.23) 

for any n-dimensional hyperplane L such that 

m^ \ProjL{ej)\ < -^-e (5.24) 

l<j<n iUU 

where Cj are the basis vectors for the x- coordinates. 

Proof of Lemma 5.2. Assume (5.14, 5.15) hold with T{yi^) replacing T{yi), where 
as before 

. logM , logMp ^ . 

lo ^ — 7-io — < ^, (5.25) 

log log 



\\{RjT{yjRj)-'\\<e<, (0 < /3 < 1) 

Mo 
10 



\{RjT{yi,)Rj)-\k,k')\ < e-^\^-^'\k,k' e J,\k-k'\ > ^ 



(5.26) 



for all 

J = [-Mo, Mo]''+^ + k, ]-M' < \k\ < M'+K (5.27) 
Recall that we are at stage z, so (H2) is satisfied. Hence 

\\T{y,)-Tiy,,)\\<Omo 

< e-"^° (5.28) 
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using (5.25). This in turn implies that (5.14, 5.15) hold. So we only need to prove 
(5.26). ((5.28) is in fact a reason for the choice of Mq in (5.5).) 

Fix / e Sjg. For type a intervals, using (H3) for yj^, the off-diagonal elements in 
the n-direction of 5" has exponential decay and \\S\\ < 0{de~'^^°). We use A. L. for 
the random Schrodinger operator eAj + Vj to obtain (5.26) as follows. 

To obtain the first estimate in (5.26), we make direct incisions in the frequency 
space. Let J be a type a interval. We require 

\±n-uj + eAj + Vj\ 
= \ ± n • u! + 

>2e-<, 

for all {n,j) G J, where /Xj are the eigenvalues of eAj + Vj restricted to the projection 
of J onto Z^. 

When n 7^ 0, this amounts to taking away a set in of measure < C»(l)e-<-Mo^+". 
When n = 0, this is satisfied for a; e X in view of (5.0). So 

||(i?jT(y,Ji?,)-i||<e<, 

by using the exponential estimates on S. This gives the first estimate in (5.26). 

To obtain the second estimate in (5.26), we use Anderson localization, i.e., on the 
probability set X defined in (4.12), for all E (here E = n ■ uj) there exists at most 

1 pairwise disjoint bad elementary regions in Z'^ of size Mq^*^ in the projection of J 
onto Z^. A resolvent series in the j direction coupled with a resolvent series in the 
n direction using the above 2 estimates and the decay property of S gives the second 
estimate in (5.26), cf. proof of Lemma 4.3. 

Hence there is a set f C Fj^ fl /, 

mes^f < e)(l)e-<M(^+^)('^+^)Mo^+^ (0 < /? < 1) 
< 

such that outside f , (5.26) hold for aU J C [-M^+i, M*+Y+'' satisfying Jn[-Mo, Mo]^x 

To prove (5.26) for type h intervals J, we use Lemma 4.1 at scale iq = [ \° ] and 

the decomposition Lemma 5.3. We illustrate this on the interval J = [— Mq, Mq]^'^'^. 
We consider the set 

S = {{to, V,e)eIxR\ (4.9) fail for T^-MoMoV+^iVio)^ i-^-, with N replaced by Mq, yi, 

replaced by y^„}, 

(5.30) 
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where / e S^^ is the same fixed interval as earher. (RecaU that x e X G \W is 
fixed.) Let T^^^iyiJ denote T'f_Mo,Mo]d+- (y^o)- Each matrix element of T^^{yi^) is a 

rational function in u, V of degree at most 2{p + 1)M'^*" {q G N"*") and linear in 9 (see 
(H4, i), (4.2-4.4, 4.7)). As before, the condition in (5.30) can be expressed in terms of 
determinants and hence polynomials in the matrix elements of T^^^y^^). This gives 

that S is semi-algcbraic of total degree at most Mq^^^'^^ M'^^o on M^'^"''-'-. 
We now localize to (a;, V) G r,;,. We consider the set 

'5' ='.5f|{(/nri„ n {DCA,c(2Mo)\nMo}) x M} C R^^', (5.31) 

where is as in (4.54) of Lemma 4.7 with N replaced by Mq, Ti^ is determined 
by the Q-equations in (2.12), which are polynomials in {uj,V) of degree at most 
2{p+ l)M'?^o(g G N+). From Lemma 4.7, (4.54), DCA,c{'^Mo)\nMo is determined 
by Mq ^^+'^'> ][j;C^{d+iy)io number of monomials of degree 1. So S' is semi-algebraic 
of total degree at most M^^+^^^o on W+^. 

By Lemma 4.1, each section S'{(jj, V) = <S'(a;(V), V) is of measure at most e~^o in 
9 {0<a< 1/2). So 

mes^+i5' < e~^o (0 < d < 1/2). (5.32) 

Our aim is to estimate ior k e , ^M^ < \k\ < M*+^ 

mes,{V| (V, V + (e + 6)^^, (V), A; • (V + (e + S)<fii, (V))) e S'}. (5.33) 

Since is a function, (/ n F^g) x M may be identified with an interval in R'^+-'^, 
say [0, 1]'^ X M, S' defined in (5.31) is a subset of (/ fl Fj^) x R, and therefore can be 
identified with a subset in [0, 1]*^ x M. For the purpose of the application of Lemma 
5.3, we identify [0, 1]" x M with [0, 1]^ x [0, l]"'^ x M and S' with S' x [0, l]"'^. Since 
T is restricted to the interval [—Mo,Mo\'^^'^, we may further restrict the interval in 
[0, 1]" X R to be 

[0, lYx[-0{Vd + ^)Mo,0{Vd + ^)Mo\ ~ [0, l]''''-W[-0{Vd + ^)Mo,0{Vdn^)Mo]. 

(5.34) 

We decompose [—0{\/d + i')Mo, 0{\/d + i')Mo\ into intervals of length 1 and iden- 
tify each of them with [0, 1]. Applying the decomposition Lemma 5.3 to each of these 
intervals and taking the union, we obtain a subset F' C Fj^ n / of measure 

mes^F' < {M^i+^^'of ■ 0{VdTl^)Mo 

< M-'/^, (5.35) 
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(Recall S' is of degree B < Mi^o(q e N+) and we take e = 200M-* in (5.22).) such 
that for all k e Z", \k\ > |M» 

mes,{V|(V, V + (e + 6)y^,,{V), k ■ {V + {e + S)^,^{V)) 

esn {((r,„\r') n / n {i)C^,e(2Mo)\OMo}) x m} 

< e-<^' (0 < (J < 1/2), (5.36) 

where Mq is as in (5.5). 

Same estimates as in (5.35, 5.36) hold when T[-Mo,Mo]''+'' replaced by T^-Mo,Mo]''+''+e^ 
ieZ'^n [-2Mo, 2Mo]'^. Therefore, there is a set T" C Tj^, n / 

mes^r" < 0{1)M-'/^ + 0{l)M^'+^^^'^+''^e-<^' 

< M-*/3 (5.37) 

(in view of the choice of Mq in (5.5) and Ca ^ 1 from (4.58)), such that outside F", 
(4.9) hold for all intervals J of the form [-Mq, Mo]'^+^ +£, £ G Z'^n [-2Mo, 2Mof, w = 
yi^ and 9 = k ■ uj,k e Z", < \k\ < M*+^ (The condition: n{DCA,c{'^Mo)\nMo} 
in (5.36) does not require additional incisions in the frequency space, as (H4, 0) holds 
starting at stage io + 1.) 

Combined with (5.29) and previous perturbation argument of replacing yi by yi^ in 
(5.28), this implies that 3r'", mes^F'" < M''/^, such that outside T'", (5.14, 5.15) 
hold for all k, <\k\< M^+i, and a fixed / G 

■C 

Letting / range over E^^, , (There can be at most (9(l)M*o of these intervals, as the 
(cd, V) -parameter space can be restricted to, say [0, l]^*^.) the total measure removed 
from is at most M*o • M~*/^ < M~*/^. Since Ti^ and Fj are at distance < < 
g-aMo from (5.3), we obtain a subset f ^ C Fj, mes^f j < M-*/^ such that (5.14, 5.15) 
hold for all k, M'/2 <\k\< M'+^ and on 

U (/n(FAf,)) (5.38) 

and hence on 

U (/n(FAf,)) (5.39) 

by (H4, iv). This proves the Lemma. □ 

Lemma 5.1 then gives that on the set in (5.39), (5.16, 5.17) hold. Clearly by per- 
turbation, (5.16, 5.17) remain valid on a M~'^*+^)'' neighborhood of (5.39), (since 
<^ g-aMo ^ ^-a{i+i)'^/^ |-jy ^YiQ choice of Mq in (5.5)), which in turn gener- 
ates a collection Ai+i of intervals in M?'^ of size 

^-(i+i) ^ such that for {u, V) G / G 
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Ai+1, (5.16, 5.17) hold. So (H4, iii) hold at stage i + 1 with a = a- M-(»+i)^ (5 > 0), 
replacing a. Moreover we have 

mes,( U (/nr,)\ U (/'nr,)) 

< mes^fi < M-*/^ (5.40) 

which will imply (H4, iv) at stage i + once we construct j/j+i and hence Fj-i-i using 
(H4, iii) at stage i + 

6. Construction of y^+i and completion of the assemblage 

Construction of yi+i 

Let N = M'+\ for (a;, V) G Ui^a^+J, define 

^i+iy = Vi+i - - {TN{yi))~^F{yi), (6.1) 

(previously (3.9)). In view of (3.1, 2.16, H4,i), this implies that Aj+ij/ is a rational 
function in (a;, V) of degree at most 



(6.2) 



0(l)Ar^+^M«*' + (2p+ 1)M«*' 
(^ e N+). 

(RecaU p < M.) So (H4, i) holds at stage i + 1. (5.16), (H4,n) give 

||A,+iy|| < M^'+^^'^Ki = 5i+^ (6.3) 

and 

\\d{A,+,y)\\ < \\dT^'\\\\F{y,)\\ + \\T^'\\\\dF{y,)\\ 

<\\T-^f\\yi\\ciKi + \\T-^\\Ki 

, .c 6.4 

where we also used (H2). Next we obtain point wise estimate on Aj+iy. From (6.1) 

|A,+iy(A;)|< ^ \T^\k,k')\\F{y{){k')\. (6.5) 

|A;'|<JV 

(2.16) gives 

\F{yi){k')\ < 0{1) Yl • • • Mk2p+i)\ 

fci+-+fc2p+i=fc' (6.6) 

< (CM)^^|/c'|('^+^)^e-"l'='l, 
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(since p < M). Substituting (5.17, 6.6) into (6.5), we then obtain 



-alfe'l 



|A,+iy(/c)| < {CMf^{ l/c'l^'^+^^^e- 

\k-k'\<iC 

+ Yl |A;'|('^+''^^e-"^l'''l+l'=-'='l^} (6.7) 

|fc-fc'|>iC 

Using (6.3) for k, such that log \k\ < i and (6.7) otherwise, we obtain 

\yi+i{k)\ < e(iogl'=l)'''-«|fe| < e-^l'^l, (6.8) 

with a = a — M~*^*+-^)^ for some ^ > independent of i, where we used the estimate 
on a just above (5.40). This shows that (H3) is essetially preserved at stage i + 1. 
Here we used the fact that a = 0{1)\ log(e + 5)\ and 0<e-Cl, 0<5-Cl. 

Since the intervals in A^+i are of size M"^*^^)'^, we may extend A^+iy to the 
entire (w, V) parameter space as follows. For any / e ^i+i, let / C / be such that 
dist(A^+i,/) ~ (l/3)M-(^+i)''. Set Ay^+i = Ay^+i on /, Ay',^, = on Af^^. Define 
a function 

1 on/ 
onA^+i. 



Define 



(6.9) 



Ay,+i = S,+iA2/^+i. (6.10) 
Ayj+i is defined on the whole (a;, V)-parameter space and satisfies 



\\dAy,+,\\ < 3M(^+i)"5,+i +M2(^+i)''^, 



(6.11) 



where the second contribution comes from (6.4). Renaming Ay^, Ayi and letting 
j/j+i = yi + Ayi, we have thus shown that (Hl-3) remain valid at stage i + 1 with a 
replacing a. 

From j/j+i, the Q-equations in (2.12) define Fj+i at most at a distance ~ 
M-^*^' < M-»/4 from F^. Clearly (5.40) implies 

mes^Fi+i n (U/eA,A ^i'€A,+, I')) < M"^/^ < M-(^+l)/^ (6.12) 
which is (H4,iv) at stage i + 
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It remains to verify the properties of in (H4,ii), stage i + From the Taylor 

series in (3.10), 

F(2/,+i) = -[(T - Tj,)[TN{y,)]-^]F{y,) + 0(l)||A,+i2/f , N = M'+\ (6.13) 

By construction, (HI), suppy^ C [-M^MY^", (2.16) gives therefore, 

suppF(2/,) C [-{2p + 1)M\ {2p + 1)MY^'' 

C [-M*+VlO,M^+VlO]^+'^ (6.14) 
= [-iV/10,iV/10]'^+^. 

So 

FiVi+i) = [Rzd+-\B{o,N)TT-^RB{o,N/io)]F{yi) + Oil^A.+iyW^, 

1 o (6.15) 

||F(y,+i)|| < ||i?z<i+.\s(o,iv)TT^'i?B(o,iv/io)||||i^(yi)|| + 0(l)||Ai+ly||^ 



where B{0,N) — [— A^, A^j^^'^, Rb{o,n)-i Rb{o,n/io) are characteristic functions. Thus 

(6.16) 



\\F{y,+,)\\<e-''/'K, + 0{l)5l, 



where we used (H2,3,4,ii). 
Similarly, 

0{l)\\dF{y,+,)\\ < \\T-m\dT\\\\Fiy,)\\ + \\dT-'\\\\F{ymT\\ 

+ \\Rz'i+'^\B{o,N)TT^^ Rb{o,n/io) II \\dF{yi) \\ 
+ \\A,+iy\\\\dA,+iy\\ 



(6.17) 



and we may take 

Ri+r = 0(1)(M2(^+i)'^«, + e-«^'^'/3«, + Si+Ji+i). (6.18) 
Summarizing (6.3, 6.4, 6.16-6.18), we have 



I Ri+i = 0(l)(M2(^+i)^/^, + e-«^^^V3^. + 5,+i5i+i). 
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(6.19) 



We start from kq, kq = 0{l){e + 5). For e + 5 small enough, (6.19) is satisfied for 
z > 1, if 



c I F, ^ lC4li I F, ^ 1/4^ + 2 (6.20) 

(H4,0) and initial input for the induction 

To ensure (H4,0) at stage i + 1, we make further incisions. (This is in order that 
Lemma 4.1 remains at our disposal at a later stage.) On r^+i, we need to eliminate uj 
such that , 

n ■ + I < e-< {n ~ O(l)Mo), 

\n ■ uJi+iWr < (0 < |n| < 2M^+i), 

where Mq is as in (5.5, H(4, 0)), = /3/0(l), 0{1) is the same expansion factor as 
in Lemma 4.6 (denoted C there), Xjj' = Hj — Hji, Hj, Hj' are the eigenvalues of the 
random Schrodinger operator eAj + Vj restricted to the myriad elementary regions 

of size M^/^^^) (the same expansion factor C(l)) in [-3Mo,3Mo]'^ (see Lemmas 4.1, 
4.2, 4.6, proof of Lemma 4.2, (4.8, 4.9), the remark after (4.9, 4.26, 4.54) and the 
definition of Qi in (H4, 0)). There are at most 0{1)M^''^ {C > 0) of such differences 
of eigenvalues. 

In view of (5.2) at stage z+ 1, the first equation in (6.21) removes a set f j+i C Fj+i, 

mes.f i+i < e-<" (0 < p" < (3') 

^ ^-[{i+if^\iogM)°'Y' (6.22) 



using (5.5) and choosing 

which is always possible. 
Since 



0(1) log M 



|a;,+i -Ui\\<5i^ M-^^'^y « ^ (6.23) 



from (5.2, 5.3, 6.20), we only need to remove a;j+i such that 

c 



|n|^ 

for M* < |n| < M'+^, which removes a set Fj-i-i C Fi_|_i, 



Oil) i±i 
mes,Fi+i < ^ « M- ^ (6.24) 
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Rename a as ccj, a as ai+i. Prom (6.8), cui+i = ccj — M"'^*+-'^)'' > Q!/2 uniformly in 
i. Combining (6.22, 6.24), we have (H4,0) at stage i + 1 and that (6.12) is preserved. 
We have thus made a complete induction from stage i to i + 1. 

7. Proof of the Theorem 

The "proof of the theorem" is now just a matter of juxtaposing sections 4, 5, 6 and 
recalling the sequence of events. We recount the spine of the argument. 

We use the modified Newton scheme in (3.9) to construct approximate solutions: 

A,+iy = y.+i -y, = -{Ti,{y,))-'F{y,), N = M'+\ (7.1) 

where T/v is T restricted to [—N^N]'^'^'^, T and F are as in (3.3-3.5, 3.1). Assume 
we have obtained the first i approximations yi, on a set of intervals M^'^ D Ai D 
• • • D Aj. To obtain we need to control {TN{yi))~^ with a further restriction to 
the new set of intervals Aj+i in (a;, V) space. This is accomplished as follows. 

To estimate TAr(yi), we cover [-M*+\ M»+i]'^+^ with the interval [-M\ M»]'^+^ 
/ and smaller intervals J = [-Mq, MqY+'' + k, /2 < \k\ < M*+^ Mq ~ (logiV)*^/^ 
as in (5.5). T^^ is "good" on Aj by using perturbation theory. The J intervals are 
divided into 2 types as in section 5, according to their distances to the Z*^ axis (see 
equations (a) and (b) between (5.13, 5.14)). TJ^ of type a is easily obtained by using 
the manifest exponential decay properties of j/j and a direct incision in the frequency 
space. The main task is to control TJ^ of type 6, which leads to further incisions in 
the frequency space, hence the new set of intervals Aj+i. 

Since \J\ <^ |/|, we may consider Tj{yig) instead of Tj{yi) for some io i as in 
(5.25). We add a parameter 9 to Tj{yi^^) and estimate the measure of the set of 9, on 
the complement of which, [Tj ]^ is "good". This is Lemma 4.1. We then use the 
decomposition Lemma 5.3 to transfer the estimate in 9 into estimates in w, giving rise 
to the new set of intervals A^+i. 

On Aj-i-i we construct yt+i according to (7.1). Using the Q-equations in (2.12), we 
obtain Pi+i. The first iq approximations are constructed by using direct e, 5 series, 
in order that Lemma 4.7 and hence Lemma 4.1 are available: zq ~ ^ log | log(e + 5)\ 

from (4.57) and the third expression in (4.9) after setting N = Nq and determining 
i, hence io- (6.20) gives the rate of convergence of this Newton scheme and hence the 
Theorem. □ 

Appendix: Localization results for random Schrodinger operators 

Random Schrodinger operator is the operator 

H ^eA + V on f{Z'^), 
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where e > is a parameter, A{i,j) = 1 if |i — j| = 1 and zero otherwise, V = {vi}i^id 
is a family of independent identically distributed (iid) random variables with common 
probability distribution g. The spectrum of H is given by 

a(H) = a(eA) + a(V) 

= [—2ed, 2ed\ + supp g, a.s. 

We summarize below the known results on Anderson localization, which are relevant 
for the present construction (cf. [DJLS, vDK, GB, GK, Mi, Si]). This is an expanded 
and more complete version of the appendix in [BW] . 

For any L G N, let ^^{i) denote any elementary region in Z'^ with diameter 2L, 
center i G Z'^ as defined in (4.10, 4.11) with Z'^ replacing U^^'^ . Let -ffA^(i) be H 
restricted to ki^{i). Let m > 0, E' G M. Ax,(i) is (m, i?)-regular (for a fixed F) if 

^(7(i?A,w) and 

|G'A.(i)(S;i,/)|<e--l^-^''l (Al) 

for all j, / G Ai;,(i), |j — / | > L/4. The following theorem is an immediate corollary 
of the corresponding theorem in [vDK] pertaining to cubes, by covering elementary 
regions with cubes and then an application of the resolvent equation (cf. Lemma B). 

Theorem A. Let I G M. be a bounded interval. Suppose that for some Lq > 0, we 
have 

Prob {for any E E I either AL^{i) or Alo(j) is {mo, E) -regular} > 1 (A2) 

for some p' > d, mo > 0, and any i,j G Z'^, |^ — j| > 2Lo 

Prob {dist {E, a{Hj^^^o))) < e~^^} < l/L^' (A3) 
for some (3 and q, < P < 1 and 

q' > Ap' + 6d (A4) 

all E with 

dist (£;,/)<^e-^', 

and all L > Lq. Then there exists a,l < a < 2, such that if we set Lk+i = [L^] + 1, 
A; = 0, 1, 2 . . . and pick m,0 < m < mo, there is Q < oo, such that if Lq > Q, we have 
that for any A; = 0, 1, 2 . . . 

Prob {for any E El I either ^Lkii) or ^Lkij) ("^) ^) regular} > 1 ^ (A5) 
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for any i,j e Z'^ with \i — j\ > 2Lfc. 
Remark. On the same probability subspace, 

dist(a(ifA,^(,)),a(i/A,j,))) > e-^S p > (A6) 

if |z — j| > 2Lk- This is part of the ingredient of the proof of Theorem A. 

Let 5 C Z"^ be an (arbitrary) finite set. Let Hs be H restricted to S. If the 
probabihty distribution is absolutely continuous with a bounded density g, we have 
the following Wegner lemma: 

Prob {dist {E, a{Hs)) < n} < Ck\S\ \\g\\oo, C > 0, k > 0. (A7) 

(A2) is verified if e is sufficiently small. (A3, 5) are provided by (A7), if ||^||oo < oo. 
More precisely, fix < /3 < 1, choose q' and hence Lq sufficiently large, then there 
exist e sufficiently small such that (A2, 3) are verified. We note from (A4) that the 
larger the q', the larger the p' could be. In view of (A5), q' can be chosen large if 
Lq is large. So p' can always be large enough by choosing e small enough to suit the 
purpose of the construction in this paper (cf. Proof of Lemma 4.1). 

Theorem A implies that for < e <^ 1, ||^||oc < oo, (t{H) has pure point spectrum 
almost surely. The pure point spectrum is dense. However it is simple [Si]. Let i/jn 
(n e Z'^) be the n*^ eigenfunction of H, then 

|V'n(j)|<C,,c.e— '1^1 (0<m'<m). 
Further improvement of technology (see [A, DJLS1,2, GB, GK]) give in fact that 

< C^P.Un,.)e-^'^'-'-'-\ (A8) 

where the centers in,w satisfy \ jn,uj\ ^ n^/*^, and is a polynomial, which only depends 
on uj. 

A resolvent estimate. 

Lemma B. Suppose A C Z*^"*"^ is an arbitrary set with the following property: for 
every x G A, there is a subset W{x) C A with x G W{x), diam {W{x)) < N and such 
that Green's function G\y^x){E) satisfies for certain t, N , A > 

\\Gwia:)iE)\\ < A (Bl) 

\Gwix){E;x,y)\<e-''' for all y e d,W {x) . (B2) 
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Here d^W{x) is the interior boundary ofW{x) relative to A given by 



d^W{x) ^ {y' eW{x)\3z e A\W{x), \z-y'\ = l}. (B3) 

Then 

\\Ga{E)\\ <2N^A 

provided AN'^e^^^ < 

See [BGS], where it is stated as Lemma 2.2, for a proof using the resolvent equation. 
See also the proof of Lemma 5.1 in section 5 of the present paper for an essentially 
identical exercise in resolvent equation. □ 
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